A Dynamic Systems Approach to Fermions and Their Relation to Spins 
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algebras to express reachable sets of pure (and mixed) states by explicit orbit manifolds. 
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I. INTRODUCTION 

Over the last decade, there has been a considerable 
experimental progress in achieving coherent control of 
ultra-cold gasses and fermionic systems p]-|fj. This is 
also of great interest in view of quantum simulation (e.g., 
0) of q uantum phase transitions pairing phenom- 

ena [10|, and in particular, for understanding phases in 
Hubbard models [Tll |. — Even earlier on, the simulation 
of fermionic systems on quantum computers had been in 
focus fl2l . Il3| . For either case, there are interesting al- 
gebraic aspects going beyond the standard textbook ap- 
proach [lJI , some of which can be found in fl5T4l7l | . Here 
we set out for a unified picture of quantum systems the- 
ory in a Lie-algebraic frame following the lines of [l8j to 
pave the way for optimal-control methods to be applied 
to fermionic systems. 

It is generally recognized that optimal control algo- 
rithms are key tools needed for further advances in ex- 
perimentally exploiting these quantum systems for sim- 
ulation as well as for computation [19 22]. In the imple- 
mentation of these algorithms it is crucial to know before- 
hand to which extent the system can be controlled. The 
usual scenario (in coherent control) is that we are given 
a drift Hamiltonian and a set of control Hamiltonians 
with tunable strengths. The achievable operations will be 
characterized by the system Lie algebra, while the reach- 
able sets of states are given by the respective pure state 
orbits. Dynamic Lie algebras and reachability questions 
have been intensively studied in the literature for qudit 
systems fhsl . [23|-[2^|. However, in the case of fermions 
these questions have to be reconsidered mainly due to 
the presence of the parity superselection rule. Hence in 
a broader sense the present work on fermions can be en- 
visaged also as a step towards quantum control theory 
for quantum simulation in the presence of superselection 
rules. 
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Apart from discussing the implications of the parity 
superselection rule in the theory of dynamic Lie alge- 
bras and of pure-state orbits, we will also treat the 
case when one imposes translation-invariance or particle- 
number conservation. Moreover, the experimentally rele- 
vant case of quasifree fermions (with and without transla- 
tion invariance) is discussed in detail. Since we interrelate 
fermonic systems with the Lie-theoretical framework of 
quantum-dynamical systems, at times we will be some- 
what more explicit and put known results into a new 
frame. 

The paper is structured as follows: In order to set a 
unified frame, we resume some basic concepts of Hamilto- 
nian controllability of qubit systems in Sec. HH since for 
comparison these concepts will be subsequently trans- 
lated to their fermionic counterparts, starting with the 
discussion of general fermionic systems in Sec. IIIII 

Then the main results are the following: In Sec. UVI we 
obtain the dynamic system algebra for general fermionic 
systems respecting the parity superselection rule (see 
Theorem @] in Subsection IIV A|) . The discussion in Sub- 
section HVB] offers an upper bound (or maximal algebra) 
in determining a system algebra by applying its symme- 
try properties (Theorem UOp . Following the same line, 
in Sec. [V] we wrap up some known results on quasifree 
fermionic systems in a general Lie-theoretic frame by 
streamlining the derivation for the respective system al- 
gebra in Proposition [T3] of Sec. |V] Corollary Q23 pro- 
vides a most general controllability condition of quasifree 
fermionic systems building on the tensor-square represen- 
tation used in (l8| . Furthermore, we develop methods for 
restricting the set of possible system algebras by analyz- 
ing their rank (i.e. the dimension of a maximal abelian 
subalgebra), see Theorem [T7l as well as Appendices [Bl and 
[Cl The structure and orbits of pure states in quasifree 
fermionic systems is analyzed in Sec. lVII leading to a com- 
plete characterization of pure-state controllability (The- 
orem [25]) ■ Sec. IVlfl is devoted to translation-invariant 
systems. For spin chains we give in Theorem 1271 the first 
full characterization of the corresponding system alge- 
bras and strengthen earlier results on short-range interac- 
tions in fl6| . The system algebras for general translation- 
invariant fermionic chains are given in Theorem [33] of 
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Sec. I VII CI The particular case of quadratic interactions 
(see Sec. IVIIDj) is settled in Theorem [Ml Corollary 1551 
considers systems which additionally observe a twisted 
reflection symmetry (or equivalently have no imaginary 
hopping terms) as discussed in [l6j . Furthermore, we pro- 
vide a complete classification of all pure quasifree state 
orbits in Sec. IVH El This leads to Sec. I VII Fl presenting 
a bound on the scaling of the gap for a class of quadratic 
Hamiltonians which are translation-invariant. Sec. IVIIII 
deals with fermonic systems conserving the number of 
particles. Their system algebras in the general case as 
well as in the quasifree case are derived in Proposition 
and Proposition |43l respectively. 

We summarize results and conclude with Sec. [IX] leav- 
ing some details and proofs to the Appendices in order 
to streamline the presentation. 



II. BASIC QUANTUM SYSTEMS THEORY OF 
iV-LEVEL SYSTEMS 

As a starting point, consider the controlled Schrodinger 
(or Liouville) equation 



[M u ,p(t)\ := -(iH u p(t) - P {t)iH u ) (1) 



driven by the Hamiltonian H u := H + X^=i u j 
and fulfilling the initial condition po : = p(0)- Here the 
drift term Hq describes the evolution of the unperturbed 
system, while the control terms {Hj} represent coherent 
manipulations from outside. Equation ([T]) defines a bilin- 
ear control system E (26| , as it is linear both in the density 
operator p(t) and in the control amplitudes Uj(t) G M. 

For a iV-level system, the natural representation as 
hermitian operators over C N relates the Hamiltonians 
as generators of unitary time evolutions to the Lie al- 
gebra u(N) of skew-hermitian operators that generate 
the unitary group U(iV) of propagators. Let L := 
{Mi, Mi, ■ ■ ■ , iH m } be a subset of Hamiltonians seen 
as Lie-algebra elements. Then the smallest subalgebra 
(with respect to the commutator [A, B] := AB — BA) 
of u(N) containig L is called the Lie closure of L writ- 
ten as (iH\, Mi, . . . , iH m )i,i e . Moreover, for any element 
M £ (Mi, . . . , M m )i\ e , there exist control amplitudes 
Uj(t) G M. with j G {!,..., m} such that 



exp(-iH) = T / exp 



-iujitfHj 



dt. 



(2) 



where T denotes time-ordering. 

Now taking the Lie closure over the system Hamilto- 
nian and all control Hamiltonians of a bilinear control 
system (E) defines the dynamic system Lie algebra (or 
system algebra for short) 



E := (M ,Mj \j = 1,2, ...,m)ue ■ 



(3) 



It is the key to characterize the differential geometry of 
a dynamic system in terms of its complete set of Hamil- 
tonian directions forming the tangent space to the time 



evolutions. For instance, the condition for full control- 
lability of bilinear systems can readily be adopted from 
classical systems |27H30j to the quantum realm such as 
to take the form of 



(M ,M j \j = l,2, 



,m)u, 



u(N) 



(4) 



saying that a iV-level quantum system is fully control- 
lable if and only if its system algebra is the full unitary 
algebra, which we will relax to su(N) in a moment. This 
notion of controllability is also intuitive (recalling that 
the exponential map is surjective for compact connected 
Lie groups), as it requires that all Hamiltonian directions 
can be generated. 

So in fully controllable closed systems, to every ini- 
tial state po the reachable set is the entire unitary orbit 
teacfifun^o) := {U p W \ U G U(7V)}. With density op- 
erators being hermitian, this means any final state p(t) 
can be reached from any initial state po a s long as both of 
them share the same spectrum of eigenvalues (including 
multiplicities). Thus the reachable set of po equals the 
isospectral set of po- 

Remark 1 . Interestingly, this notion is stronger than the 
requirement that from any given (normalized) pure state 
one can reach any other (normalized) pure state, since it 
is well known (23l - t25| that for N being even, all rank-one 
projectors are already on the unitary symplectic orbit 



teach(|Vo>(^o|) = {K\i> )(ip \Kl | K G Sp(7V/2)} 
= {U\4, )(ip \tf\U€SXl(N)} 



(5) 



and Sp(iV/2) is a proper subgroup of SU(7V). 

In general, the reachable set to an initial state po of a 
dynamic system (E) with system algebra jjs is given by 
the orbit of the dynamic (sub)group := cxp(gs) G 
U(N) generated by the system algebra 

reacf) E (po) := {Gp G t | G G G s = exp( 0E )}. (6) 

Thus the system algebra qy. can be envisaged as the 
fingerprint encoding all the dynamic properties of a dy- 
nami c sy stem E. Via the respective reachable sets (see, 
e.g., [18j]) it is easy to see that a coherently controlled dy- 
namic system E^ can simulate the dynamics of another 
system E# if and only if the system algebra qy, a of the 
simulating system E^ encompasses the system algebra 
of the simulated system Eg, 



(7) 



In [18| , we have analyzed the possibility of quantum sim- 
ulation with respect to the dynamic degrees of freedom 
and have given a number of illustrating worked examples. 

Next we describe dynamic symmetries of bilinear con- 
trol systems whose Hamiltonians are given by m :— 
{iH„} = {Ma, Mi, . . . , M m }. The symmetry operators 
s are collected in the centralizer of m in u(iV): 

cent(m) :={se u{N) \ [s,M] = Vii? G m}. (8) 
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More generally, let S' denote the commutant of a set S 
of matrices in the set of all complex matrices. By Ja- 
cobi's identity [[a, b] , c] + [[b, c] , a] + [[c, a] , b] =0 one gets 
two properties of the centralizer pertinent for our con- 
text: First, an element s that commutes with all Hamil- 
tonians a, b G m also commutes with their Lie closure 
S := (m) Lic (i.e. cent(m) = cent(g s )), as [s,a] = and 
[s, b] = imply [s, [a, b]} = 0. Second, for any u G u(N), 
[s\,u] = and [s2,u] — imply [[si,S2],w] = 0, so the 
centralizer forms itself a Lie subalgebra to u(N) consist- 
ing of all symmetry operators. By construction, cent(fjs) 
is also a normal subalgebra or an ideal of u(N) observing 
[cent(fl S ),u(JV)] C cent( 0s ). 

Likewise one can describe the symmetries to a given 
set py, of states by its centralizer 

centos) := {s G u(N) \ [s, p] = Vp G p s } (9) 
= cent((ps)a), 

where (-)r denotes the real span. Clearly, cent(ps) C 
u(iV) generates the stabilizer group to the state space ps 
of the control system (£). 

Since in the absence of other symmetries the identity 
is the only and trivial symmetry of both any state space 
Py, as well as any set of Hamiltonians and their respec- 
tive system algebra qy, one has cent(gs) = centos) = 
{i XIn I A G R} =: u(l). So there is always a trivial sta- 
bilizer group U(l) := {e^lw |0 G R}. This explains 
why the time evolutions generated by two Hamiltonians 
Hi and H2 coincide for the set of all density operators if 
(and without other symmetries only if) Hi —H2 = Al. As 
is well known, by the same argument, in time evolutions 

p(t):=U(t)ptf(t) = Ad uit) (p ) (10) 

following Eq. (fTJ) , one may take U(t) := exp(-itH) 
equally well from U(iV) or SU(A). Thus henceforth we 
will only consider special unitaries (of determinant +1) 
generated by traceless Hamiltonians iH v G su(N), since 
for any Hamiltonian H there exists an equivalent unique 
traceless Hamiltonian H := H — tr(£T)ljv generating 
a time evolution coinciding with the one of H [3l| . 

However, the above simple arguments are in fact much 
stronger, e.g., one readily gets the following statement: 

Lemma 2. Consider a bilinear control system with sys- 
tem algebra Qy on a state space py- Let iHi G 0s an d 
iH2 G u(N) while assuming that [Hi, (py)r\ Q i(pY,)is. for 
all iHi G Qy, *-e., operations generated by Qy map the set 
(Ps)r into itself. Then the condition 

e -iH ltpe i Hl t = e -i{H 1+H2) t pe i{H l+ H 2 )t vt G R)p g ps 

(11) 

is equivalent to iH2 G cent(ps)- 

Proof. Using the formula e tA Be~ tA = cxp[adtA(B)] = 
J2kLo t k /k\ ad^(-B) we show that Eq. pT|) is equivalent 
to condition (a): ad k Hi {p) = ad^ 1+J? (p) for all non- 
negative integer k and all p G (ps)r- Moreover, (a) 



implies condition (b): (ad# 2 o&d Hl )(p) = for all non- 
negative integer k and all p G (/Os)r, as [Hi, ad^ 1 (p)] = 
[Hi + H 2 M k ^\ H2 {p)} = [Hi +H 2 ,ad k H - 1 (p)}. Also, (a) 
follows from (b) due to ad k Hl (p) = [Hi+ H 2 ,ad h I f 1 1 (p)} = 
[Hi + H 2 , [Hi + H 2 ,ad k H - 2 (p)}} =■■■ = ad k Hi+H2 (p). Ap- 
plying [Hi, (ps)r]Ci(pe) k to (b) completes the proof. □ 

Therefore, let us consider a pair of Hamiltonians 
iHi , iHj, G Qy, (fulfilling the conditions of Lemma [5]) 
as equivalent on the state space p-%, if the difference 
ii?2 := i(Hi — H3) between them falls into the centralizer 
cent(p E ). 

III. FERMIONIC QUANTUM SYSTEMS 

In this section, we fix our notation by recalling basic 
notions for fermionic systems. In the first subsection, 
we discuss the Fock space and different operators acting 
on it as given by the creation and annihilation operators 
as well as the Majorana operators. We point out how 
the Lie algebra u(2 d ) of skew-hermitian matrices can be 
embedded as a real subspace in the set of the complex 
operators acting on the Fock space. In the second sub- 
section, we focus on the parity superselection rule and 
how it structures a fermionic system. 

A. The Fock Space and Majorana Monomials 

The complex Hilbert space of a d-mode fermionic sys- 
tem with one-particle subspace C d is the Fock space 

= c e c d © (c d a c d ) © • • • © ^/\ c^J . 

Given an orthonormal basis {e,}|_j of C d , the Fock vac- 
uum £1 := 1 (= 1 © © • • • © 0) and the vectors of the 
form A ei 2 A • • • A e^ (with ii < i 2 < ■ ■ ■ < and 
1 < k < d) form an orthonormal basis of J-"(C d ). Note 
that F(C d ) is a 2 d dimensional Hilbert space isomorphic 

to ©f =1 C 2 (S C 2 "). 

The fermionic creation and annihilation operators, /J 
and f p act on the Fock space in the following way: f^Cl = 
e p , f p Q = 0, p p e q — e p A e q , and f p e q = S pq ; while in 
the general case of 1 < £ < d, their action is given by 
/t(e gi A e q . 2 A • ■ ■ A e qe ) = (e p A e qi A e q2 A • • ■ A e qt ) and 

f p (e qi /\e q2 A- ■■Ae qe ) = ELiC -1 )* %ih e n A " ' - Ae <i(k-i) A 
e 9(fc+i) A • • • A e 9f . By their definition, these operators 
satisfy the fermionic canonical anticommutation relations 

{fl f\) = {f P , f q ) = and {/t, /,} = S pq l, 

where {A, B} := AB + BA denotes the anticommutator. 
Moreover, every linear operator acting on T(C d ) can be 
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be written as a complex polynomial in the creation and 
annihilation operators. 

Another set of polynomial generators acting on the 
Fock space is given by the 2d hermitian Majorana op- 
erators m 2p -i := f p + P p and m 2p := i(f p - /+), which 
satisfy the relations (k, I <G {1, . . . , 2d}) 



{m kl m e } = 26 u l. 



A product m q 



of k > Majorana operators 



is called a Majorana monomial. The ordered Majorana 
monomials with q\ < q 2 . . . < q k form a linearly inde- 
pendent basis of the complex operators acting on P(C rf ). 
Each Majorana monomial acting on d-mode fermionic 
system can be identified with a complex operator acting 
on a chain of d qubits via the Jordan- Wigner transfor- 
mation (33435| which is induced by 



m 2p -i 



H> Zi 



z<e>x® ii 



1 



p-i 



d—p 



and 



m 2 p <— > Z 1 



where the following notation for the Pauli matrices X := 

(? o)> Y:= (? ~o)> and z := (0-1) is used - 

Now we highlight the real subspace contained in the 
set of complex operators acting on the Fock space .F(C d ) 
which consists of all skew-hermitian operators and which 
forms the real Lie algebra u(2 d ) closed under the commu- 
tator [A, B] := AB — BA and real-linear combinations. 
More precisely, u(2 d ) is generated by all operators 



L(M) := ~w{M)M 



(12) 



where M denotes any ordered Majorana monomial and 



w(M) :-- 



i if [deg(M)mod 8] € {0,1}, 

1 if [deg(M)mod 8] G {2,3}, 

-i if [deg(M ) mod 8] G {4,5}, 

-1 if [deg(M)mod 8] G {6,7}. 



(13) 



Similarly, one obtains a basis of su(2 d ) by excluding — 11. 



a commuting observable in turn implies that a super- 
position of pure states from different blocks of a block- 
diagonal decomposition w.r.t. the eigenspaces of this ob- 
servable are equivalent to an incoherent classical mixture. 

The parity super selection rule identifies among the op- 
erators acting on F(C d ) the physical observables Mp as 
those that do commute with the parity operator 



2d 



P 



(14) 



fe=i 



where the adjoint action of P on a Majorana monomial 
isgivenasPm fel TOfe 2 • • ■m kl P^ 1 = (-l) e m kl m k2 ■ ■ -m kr 
These physical operators are also exactly the ones that 
can be written as a sum of products of an even number of 
Majorana operators (as P contains all Majorana opera- 
tors whereof there exist an even number), and are there- 
fore also called even operators for short. If the parity is 
the only non-trivial symmetry, we obtain H^. = (1,P), 
where the bracket stands for the complex-linear span. 

Now we will discuss why the set of all physical 
fermionic states Pf consists similarly of all density op- 
erators that commute with P, notably p' F = (1,P). As 
we will show, the parity superselection rule induces a de- 
composition into a direct sum of two irreducible state- 
space components exploiting M' F D p' F = (1,P). Re- 
call that P 2 = 1 and the eigenspaces to the eigenval- 
ues +1 and —1 are indeed of equal dimension, as there 
are exactly 2 2d ~ 1 even operators which map the vac- 
uum state fi into the +1 eigenspace of P. Note that 



Pe qi A e Q2 A 



A e qe = (-l)'e 



<?1 ^ ^*?2 ^ 



A e qe . Thus 



the Fock space can be split up as a direct sum of two 
equal-dimensional eigenspaces of P, called the positive 
and negative parity subspaces (for clarity observe (13|): 



J-(C d ) 



e a- 



e 



A 1 



xld 



Now we may write P 2 = 1 = P+ + P- with the 
orthogonal projections P + := 5(1 + P) and P_ := 
i(l — P) projecting onto the respective subspaces. Any 
physical observable (i.e. even operator) A has a block- 
diagonal structure with respect to the above splitting, 
i.e. A = P+AP+ + P_AP-. This follows, as the re- 
quirement [A,P] = |[A,P+] = -\[A,PJ\ = en- 
forces P + AP- = P_AP + = for any operator A = 
P+AP+ + P+AP- + P-AP+ + P_AP-. We obtain 



B. Parity Superselection Rule 



Tr(pA) = Tr{pP + AP + 
= Tr[(P+pP+ 



f pP_AP_) 
P-pP-)A\. 



(15) 



An additional fundamental ingredient in describing 
fermionic systems is the parity superselection rule. Su- 
perselection rules were originally introduced by Wick, 
Wightman, and Wigner [36j (see also [13, [38[)- These 
rules, in the finite-dimensional definition of Piron (39|, de- 
scribe the existence of non-trivial observables that com- 
mute with all physical observables. The existence of such 



Hence physical observables cannot distinguish between 
the density operator p and its block-diagonal projection 
to P + pP + + P-pP^ (which is always an even density 
operator). In this sense, a physical linear combination (a 
formal superposition) of pure states from the positive and 
negative parity subspaces is equivalent to an incoherent 
classical mixture. Equation (| 1 5[) also shows that without 



5 



loss of generality we can restrict ourselves to even density 
operators and regard only those as physical. 

Finally, we would like to recall three further aspects of 
the parity superselection rule. First, without the parity 
superselection rule, two noncommuting observables act- 
ing on two different and spatially-separated regions would 
exist which would allow for a violation of locality (e.g., 
by instantaneous signaling between the regions) . Second, 
the parity superselection rule, of course, does not apply if 
one uses a spin system to simulate a fermionic system via 
the Jordan- Wigner transformation. This system respects 
locality, since the Majorana operators vrik are — in this 
case — localized on the first [(fc+1) div2] spins; two non- 
commuting Majorana operators are therefore not acting 
on spatially-separated regions. Third, the parity supers- 
election rule also affects the concept of entan glem ent as 
has been pointed out and studied in detail in 



glement 



IV. FULLY CONTROLLABLE FERMIONIC 
SYSTEMS 



Here we derive a general controllability result for 
fermions obeying the parity superselection rule. We es- 
tablish a double-centralizer theorem for the dynamic sys- 
tem algebras (not to be confounded with von Neumann's 
double-commutant theorem) and conclude by an example 
illustrating full controllability for a fermionic system with 
quadratic interactions and a single fourth-order term. 



A. System Algebra 

In the case of qubit systems mentioned in Sec. UU two 
Hamiltonians generate equivalent time evolutions if and 
only if they differ by a multiple of the identity. This 
condition can readily be modified for the fermionic case 
such as to match the parity-superselection rule as well. 

Corollary 3. Let H\ and H2 be two physical fermionic 
Hamiltonians, i.e., even hermitian operators acting on 
J-(C d ). Then by Lemma\^the equality 



and moreover, H and H are equivalent and generate 
the same time evolution. If necessary, we can restrict 
ourselves to the set of Hamiltonians satisfying Eq. (JTTj) . 
These elements decompose as H = H + H _ , where H + 
and if_ are generic traceless hermitian operators each 
acting on a 2 rf_ ^dimensional Hilbert space. We explic- 
itly define the linear space F^ of physical fermionic Hamil- 
tonians as generated by the basis of all even Majorana 
monomials without the operators 1 and P, ensuring that 
is traceless both on H + and H_ . — We summarize our 
exposition on fully controllable fermionic systems in the 
following result: 

Theorem 4. The Lie algebra corresponding to the phys- 
ical fermionic (and hermitian) Hamiltonians F d is 



C d := su^ ) ® su( 2 ) 



(18) 



The most general set of unitary transformations gener- 
ated by Cd is given as the block- diagonal decomposition 
SU(2 d ' 1 )©SU(2 d - 1 ). Henceaset{H ,H u H 2 ,...,H m } 
of hermitian Hamiltonians defines a fully controllable 
fermionic system if and only if 



-iHit iH x t 



e -iH 2 t p( ,iH 2 t 



(iH ,iH 1 ,...,iH m ) Uc =su(2 d - L )®su(2 cl - 1 ) . (19) 

Remark 5. For Lie algebras, £i + £2 will denote only an 
abstract direct sum without refering to any concrete re- 
alization. We reserve the notation ti © t-i to specify a 
direct sum of Lie algebras which is (up to a change of 
basis) represented in a block-diagonal form ( {l {2 ) . 

Proof. It follows from Sec. IIIII that F^ commutes with 
P and that the matrix representation of F^ splits into 
two blocks of dimension 2 d ~ 1 corresponding to the + 
and — eigenspaces of P. As the center of F^ is given 
by W' d H Fd = (1,P) n F d = {0}, the Lie algebra F d 
is semisimple. As there are exactly 2 2<i ~ 1 — 2 linear- 
independent operators in F^, the system algebra could 
be su(2 d ~ 1 ) ffisu(2 d ~ 1 ). And indeed, all other system al- 
gebras are ruled out as the subalgebras acting on each of 
the two matrix blocks would have a smaller Lie-algebra 
dimension than su(2 d_1 ). □ 



B. Examples and Discussion 



holds for all even (physical) density operators pf with 
p' F — (1,P) in the sense that H\ and H2 generate the 
same time- evolution, if and only if H 2 — H\ = Xl + fxP = 
(A + n)P+ + (A - p)P- with A, fi e M. 

This also implies that for any physical fermionic Hamil- 
tonian H , there exists a unique Hamiltonian 



H := H 



tr(P + HP + 
dim 



P + - 



tr(P-gP-) 
dim P 



(16) 



that is traceless on both the positive and the negative 
parity subspaces, i.e., 



tr(P + HP + ) = tr(P_PP_) = . 



(17) 



We start by providing an example realizing a fully con- 
trollable fermionic system: 

Proposition 6. Consider a fermionic quantum system 
with d > 2 modes. The system algebra Cd = su(2 d ^ 1 ) © 
su(2 d_1 ) of a fully controllable fermionic system can be 
generated using the operators W\ := L(v\), W2 '■= L(v>2), 
W3 := L(v%), and u>4 := L(vi) with the map L as defined 
in Eqs. (|12[) and (|13p . where 



(20a) 



d-l 

vi := -m 2p -im 2 p + 2 + rn 2 pm 2 p + i, 
P =i 

v 2 := mxn%2, W3 := m 2 m^, va '■= m\m2 , mz'm^. (20b) 
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Proof. All the generators are contained in Cd by ar- 
guments of Subsection IIV A[ One can verify with ex- 
plicit computations that the given generators yield any 
Majorana operator of even degree besides L(riq=i m <?)- 
Counting the operators generated, we obtain that the 
dimension of the subalgebra generated is equal to the 
dimension of Cd- Therefore, the proposition follows. □ 

The proof also implies that all the operators generated 
commute with Yll=i m q = P/i d [cf- Eq. (|14p ] (and the 
identity operator 1). In addition, all operators commut- 
ing simultaneously with all elements of Cd can be written 
as a complex-linear combination of 1 and P. We thus 
obtain a partial characterization of full controllability in 
fermionic systems: 

Lemma 7. Consider a fermionic quantum system with 
d > 2 modes. A necessary condition for full controllabil- 
ity of a given set of hermitian Hamiltonians H v is that 

{iH v y = (i,p). 

One can expect that the condition of Lemma [7] is not 
sufficient under any reasonable assumption by applying 
counterexamples from spin systems in [l8j . These coun- 
terexamples could be lifted to fermionic systems by pro- 
viding the explicit form of the embeddings from su(2 d_1 ) 
to the first and second component of the direct sum 
C d =5u(2 d " x ) ©su(2 d - 1 ). 

We guide the discussion in a different direction by em- 
phasizing that the property {iH v }' = (1,P) does not 
determine the system algebra uniquely. We define the 
centralizer of a set B C su(fc) in su(fc) (e.g. k = 2 d ) as 

cent 5u(fc) (5) := {g e su(fc) | [g, b] = for all b e B}. 

We consider the algebras Cd = su(2 d_1 ) ©su(2 d_1 ) and 
s[u(2 d_1 ) © u(2 d_1 )], where the latter algebra is isomor- 
phic to su(2 d_1 ) + su(2 rf_1 ) + u(l) and contains the 
additional (non-physical) generator L(n<j=i m q)- Note 
that ceni su(fe )(£d) = cent 5u(fc) (s[u(2 rf - 1 ) © u(2 rf - 1 )]) = 
L(IIg=i m q)i i- e -j the centralizers of both algebras are 
equal. However cent su(fc ) [Ldlgfa m g)] = s[u(2 d_1 ) © 
u(2 d " 1 )] ^ su(2 d - 1 ) ffisu(2 d " 1 ). In particular, we have 
C d 7^ cent 5U ( fe )(cent 5u(fc )(£ d )), and C d does not fulfill the 
double-centralizer property. A more general incarnation 
of this effect is the following formal consideration: 

Proposition 8. Let g denote a subalgebra of su(fc). 
There exists a set A C 5u(k) such that g = centra,) (A), 
if and only if cent 5u(fc) (cent su(fc) (g)) = g. 

Proof. First, let us assume the existence of the set A. As 
cent su(fe )(cent 5u(fc) (cent 5u(fe) (A))) = cent 3U ( fc )(A]holds for 
any set A, which can also be inferred from [4J, Propo- 
sition 6.1.3.1(iii)], we obtain cent 6U ( fe ) (cent su ( fc ) (gj) = g. 
Second, we assume that cent 5u ( fc ) (cent su ( fc ) (g)) = g holds. 
We choose A := cent^M (g) and verify its existence. □ 



To further analyze the influence of symmetry proper- 
ties on the system algebra, we recall some elementary 
representation theory (see, e.g., Theorem 1.5 of (iij): 

Proposition 9. Consider a completely reducible complex 
matrix representation $(;?) of a group G, where k is the 
degree of $. Let comm( ( I > ) = <I > ' denote the commutant 
algebra of all complex k x k-matrices simultaneously com- 
muting with $((?) for g € G. Then, <£(<?) is equivalent to 

w 

[l e . , 

where <pj denote for j G {1, . . . , w} distinct inequivalent 
irreducible complex matrix representations of G with de- 
gree kj , occuring with multiplicity ej in <!> . In particular, 

(a) dimcomm($) = J2j=i e % 

(b) dimcenter(comm($)) = w, 

(c) k = J2J=1 k j e j- 

Obviously, the same is true for representations of a 
compact Lie group or its Lie algebra. Given a subalgebra 
g of su(k) (or respectively of u(fc)) and a representation 
$ of g with degree k, we discuss the easiest case of Propo- 
sition [5] where w = 1 and ei = 1. Hence, <!> is irreducible 
and g is an irreducible subalgebra of su(fc) (or respec- 
tively of u(k)). But g is not necessarily equal to su(fc) (or 
respectively to u(fc)). Irreducible simple subalgebras of 
su(fc) were studied extensively in this regard in Ref. [l8[ . 
Note that the irreducible subalgebras of u(k) are of the 
form g or g+u(l) where g denotes any irreducible subalge- 
bra of su(fc) (cf. pp. 27-28 and p. 321 of (H). — A slight 
generalization is given by the case of an abelian commu- 
tant algebra, i.e. dim comm(<l>) = dim center(comm($)) 
and ej = 1 for all j £ {1, . . . , w \. O ne may thus apply 
the spectral theorem (see, e.g., |4q - |48l |) simultaneously 
to all the elements of the commutant algebra: 

Theorem 10. Consider a Lie algebra g C su(fc) and 
its representation $ of degree k. Assume that the cor- 
responding commutant algebra C = comm( < I ) ) is abelian. 
One obtains that g is a subalgebra o/s[©^™ u(kj)) and 

it is equivalent to s[ffij™ c 0j], where k = X)j™ C k j an d 
gj are irreducible subalgebras ofu(kj). Furthermore, one 
finds kj = dim(fj), where Pj are the orthogonal projec- 
tion operators given by the joint spectral decomposition 

of C with Y^j=l C Pj = ■"■*! an d PiPj = for i ^ j. If 
g is the maximal Lie algebra with these properties, then 
g = 5 [(Bf™ 1 c u(k J )}. 

Using Proposition |H] one can directly characterize a 
maximal Lie algebra g contained in su(fc) which is defined 
by all its symmetries including cases where the commu- 
tant to g is not necessarily abelian. Observe the notation 
of Remark 

Theorem 11. Consider a Lie algebra g C su(A:) and its 
representation <f> of degree k. Let C = comm( ( E > ) denote 
the commutant of g. If g is the maximal Lie algebra 
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with these properties, then g = s[ 
dim[centet(C)] and Y^j=i kj — 



i=i u (kj)] where lo 



Proof. Using Proposition[n](and its notation) one obtains 
that g is equivalent to (Bj = i [l ej <£> </>j(g)\ ■ Therefore, g 
is a subalgebra of s[+^ =1 u(/cj)] with Y^j=i kj — The 
maximality of g completes the proof. 

In a dual approach, one could start from a set S of 
symmetries of g. Due to the maximality of g, the set S 
has to comprise all symmetries of q. Next, one can apply 
Proposition[5]to the subalgebra of su(k) generated by the 
linear span intersected with su(fc), i.e. (S) n su(fc). The 
theorem then follows directly using Schur's lemma and 
the maximality of g. □ 

The reader familiar with the double-commutant theo- 
rem in algebraic quantum mechanics will wonder about 
the different power of symmetries for characterizing al- 
gebras of observables on the one hand and Lie algebras 
on the other: a von-Neumann algebra A is entirely de- 
termined by its commutant A', since A" = A |49l. [50|. 
In this sense, there is a duality between the algebra A 
and its commutant A' encapsulating all symmetries. On 
the other hand, consider the illustrative case of an ir- 
reducible Lie subalgebra g of su(k) [5l|, where the cen- 
tralizer cent su (fc) (g) is trivial, i.e. zero. This centralizer 
is shared with all irreducible Lie subalgebras of su(fc). 
So in turn, the double centralizer in su(k) to all these 
subalgebras is su(k) itself. We thus obtain the following 
corollary to Proposition [5] and Theorem [TU1 where the 
double centralizer gives a maximality criterion ensuring 
that an irreducible g C su(fc) is in fact fulfilling g = su(k) 

Corollary 12. Let g denote an irreducible subalgebra 
of su(k), i.e. cent su ( fc )(g) = {0}. Then one finds that 
cent su(fe )(ccni«, u(fc) (g)) = g if and only ifg = su(k) JM]- 

Summarizing the general case, the symmetry proper- 
ties of a Lie algebra g C su(fc), as given by its commutant 
w.r.t. a representation of g, do not determine the Lie al- 
gebra g uniquely. Yet the commutant allows us to infer 
a unique maximal Lie algebra contained in su(fc), which 
is (up to an identity matrix) equal to the double com- 
mutant of g, but in general not to g itself. Although all 
representations of compact Lie algebras, such as su(k) 
and its subalgebras, are completely reducible, the situa- 
tion for Lie algebras also differs from the case of associa- 
tive algebras: here complete reducibility of a representa- 
tion implies the double-commutant theorem (see Theo- 
rem (3.5.D) of [54j or Theorem 4.1.13 of [HBj). whereas 
the double-commutant theorem does not apply to Lie al- 
gebras as discussed above. 



V. QUASIFREE FERMIONS 

Here we present the dynamic system algebras for 
fermions with quadratic Hamiltonians. For illustration, 



the relation to spin chains with Heisenberg-XX interac- 
tion is given, which allows to adapt our tensor-square 
criterion for full controllability from spin systems [l8| to 
quasifree fermionic systems. 



A. Quadratic Hamiltonians 

A general quadratic Hamiltonian of a fermionic system 
can be written as (cf. jlH |56l - (59| ) 

d 11 

H = A pq (f p : f q — 5p q ^) + -Bpqf^f^ — -B pq f p fq, 

p,q=l 

(21) 

where the coupling coefficients A pq and B pq are complex 
entries of the d x d-matrices A and B. respectively. The 
canonical anticommutation relations and the hermicity of 
H require that A is hermitian and B is (complex) skew- 
symmetric. The terms corresponding to the non-zero ma- 
trix entries of A and B are usually referred to as hopping 
and pairing terms, respectively. Related parametriza- 
tions for quadratic Hamiltonians are discussed in Ap- 
pendix |XJ 

In the Majorana monomial basis, the quadratic Hamil- 
tonian H can be rewritten such that 



2(/ 



with 

T=i 

2 



iH 



i 

--m k rne 



(22) 



-°i ;)+*«(*)*(; ^ 



-3m(A) ® ( X + 3m(B) ® ( ^ Jj 



The properties of A and B directly imply that the matrix 
T is real and skew-symmetric. Using the formula 

[m p m q ,m r m s ] = -A(S ps S qr l - S qs 5 pr l) 
+2(5 ps m q m r —Sp r m q rn s +S qr rnprn s —5q S rnprn r ) (23a) 

= 5p S (m q m r — m r m q ) — 5 pr (m q m s — m s m q ) 
+6 qr (m p m s — m s m p ) — S qs (mpm r — m r m p ) (23b) 

one can easily verify that the space of quadratic Hamil- 
tonians is closed under the commutator. To sum up, we 
have established the well-known Lie homorphism from 
the system algebra generated by a set of quadratic Hamil- 
tonians (whose control functions are given by the matrix 
entries of A and B) onto the system algebra so (2d) rep- 
resented by the entries of T (cf. pp. 183-184 of |35|): 

Proposition 13. The maximal system algebra for a sys- 
tem of quasifree fermions with d modes is given by so(2d). 

Proof. Let the map h transform the Majorana mono- 
mial — ^(m p m q —m q mp) into the skew-symmetric ma- 



trix e r , 



where e pq has the matrix entries [e pq ] z 
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SpuSqy We show that h is a Lie-homomorphism assuming 
p =/= q and r =/= s in the following, while the case of p = q 
or r = s holds trivially. Note that \{m p m, q — m q m p ) = 
m p m q . It follows from Eq. (|23b[) that h([— ^(m p rriq — 
m q m p ),-\(m r m s -m a m r )}) = [{e pq -e gp ), (e rs -e sr )] = 
[h(—^(m p m q — m q m p )), h(—^(m r m s — m s m r ))]. □ 

B. Examples and Explicit Realizations 

We start by showing that the full system algebra 
so(2d) of quasifree fermions can be generated using only 
three quadratic operators, namely w\ = L(v\), w 2 = 
L(v 2 ), and u>3 = L(v%) from Eq. (f20|) where v\ = 

S p =i ~ m 2p-im2 P+2 + m 2p m 2p+ i, v 2 = mim,2, and v 3 = 
m 2 m 3 . The Jordan- Wigner transformation maps these 
generators respectively to -f Ep=i Vp+i + Y p Y p+ i, 
— |Zi, and — |XiX2, where operators as (e.g.) Zi are 
defined as Z ® I ® • • • (g) I- 

Lemma 14. Consider a fermionic quantum system with 
d > 2 modes. The system algebras t\ and t 2 generated by 
the set of Lie generators {wi,w 2 } and {w\, W2, 103} con- 
tain the elements L{a p ) with a p := m 2p -\m 2p for all p £ 
{1, . . . , d} as well as L(b p ) with b p := ~m 2p -im 2p+2 + 
m 2p m 2p+ i and L(c p ) withc p := m 2p -im 2p+ i+m 2p m 2p+2 
for all p £ {1, . . . , d — 1}. 

Note that the elements L(a p ), L(b p ), and L(c p ) are 
mapped by the Jordan- Wigner transformation to the 
spin operators —iZ p /2, —i(X p X p+ i + Y p Y p+ i)/2, and 
-i(X p Y p+1 - Y p X p+1 )/2, respectively. 

Proof. We compute the commutators W4 := —L(c,\) — 
[v)2,Wi], w 5 := L(bi) = [104,102], and w 6 := L(a 2 ) = 
[105,104] — w 2 . We can now reduce the problem from 
d to d — 1 by substracting 105 from w\. The cases of 
d £ {2,3,4} can be verified directly and the proof is 
completed by induction. □ 

This proof also yields an explicit realization for the 
algebra so (2d) while providing a more direct line of rea- 
soning as compared to our proof of Theorem 32 in [l8[ . 

Theorem 15. Consider a fermionic quantum system 
with d > 2 modes. The system Lie algebra t 2 generated 
by {ioi, 102, 1^3} is given by so (2d). 

Proof. The cases of d £ {2, 3, 4} can be verified directly. 
We build on Lemma [T4l and remark that % 2 Q so (2d) as 
it is generated only by quadratic operators (see Propo- 
sition 1 13|) . We compute in the Jordan- Wigner picture 
w 7 := -i(YiY 2 - Y 2 Y 3 )/2 = [103, [w 3 , wi]], and w 8 ■= 
—1X2X3/2 = L(b 2 ) — (105 — W3 — 107). This shows by 
induction that so(2d) D t 2 2 u(l) + so(2d - 2). As 
u(l) + so(2d — 2) is a maximal subalgebra of so(2d) (see, 
e.g., p. 219 of [60]), one obtains that l 2 = so(2d). Al- 
ternatively, one can explicitly show that t 2 consists of 
all quadratic Majorana operators, which combined with 
Proposition Q2] would also complete the proof. □ 



Note that the generators wi, 102, and W3 can be de- 
scribed using the Hamiltonian of Eq. (|21[) while keeping 
the control functions given by the matrix entries A pq and 
B pq in the real range, see Appendix |A1 for details. This 
also provides a simplified approach to Theorem 32 in [l8| , 
where only the real case was considered: 

Corollary 16 (see Theorem 32 in [H[)- Consider a 
control system given by the Hamiltonian components of 
Eq. pip . The control functions are specified by the ma- 
trix entries A pq and B pq which are assumed to be real. 
The resulting system algebra is so (2d). 

Relations of quasifree fermions to spin systems will be 
analyzed in Appendix [Bj — Next we treat the case of the 
algebra u(d). 

Theorem 17 (see Lem. 36 in [18]). Consider a fermionic 
quantum system with d > 2 modes. The system Lie alge- 
bra 61 generated by {wi,w 2 } is given by u(d). 

Here we just sketch ideas for the proof of Theorem [T7] 
while leaving the full details to Appendix [C] Our meth- 
ods exploit the detailed structure of the appearing Ma- 
jorana operators while being more explicit than in (l8| 
and avoiding obstacles of the spin picture. Building on 
the notation of Lemma [TH we show that the elements 
L(a p ) with 1 < p < d together with the elements L(b p ^) 
with b£> := -m 2p _im 2p+2l + m2 P m2 P +2i-i and L(cp ) ) 

with c p := m 2 p_iTO2p+2i-i + m 2p m 2p+2 i where p, i > 1 
and p + i < d form a basis of i\. One obtains that 
dim(fii) = d + (d — l)d = d 2 . Furthermore, the elements 
L(a p ) form a maximal abelian subalgebra and the rank 
of £i is equal to d. We limit the possible cases further by 
showing that 61 is a direct sum of a simple and a one- 
dimensional Lie algebra. A complete enumeration relying 
on tables in [6(j completes the proof. 

C. Tensor-Square Criterion 

Consider a control system of quasifree fermions which 
is represented by matrices T„ in the form of Eq. <\22\ . For 
more than two modes (i.e. d > 3), we can efficiently de- 
cide if the system algebra is equal to so (2d). Recall that 
the alternating square Alt 2 (</>) and the symmetric square 
Sym 2 (0) of a representation 4> are defined as restrictions 
to the alternating and symmetric subspace of the tensor 
square (j)® 2 = ® l dimW + l d i m (0) ® 

Theorem 18. Assume that t is a subalgebra of so (2d) 
with d > 3 and denote by $ the standard representation 
o/so(2d). Then, the following are equivalent: 

(1) t = so(2d). 

(2) The restriction of Alt 2 <£> to the subalgebra t is irre- 
ducible and the restriction o/Sym 2 $ to 4 splits into two 
irreducible components. Each irreducible component oc- 
curs only once. 

(3) The commutant of all complex matrices commuting 
with the tensor square ($|{)® 2 oft has dimension three. 
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Proof. Assuming (1), condition (2) follows from the for- 
mulas for the alternating and symmetric square of so (2d) 
with d > 3 given in its standard representation 0(1,0,. •• ,o) 
[where (1,0, ••• ,0) denotes the corresponding highest 
weight]: The alternating square is given as Alt 0(i,q,o) = 

0(o,i-i) for so ( 6 ) and Alt2 0(i,o.o,-,o) = 0(o,ul--,o) for 
so (2d) in the case of d > 3 (cf. Table 6 in [61| or Ta- 
ble X in [Hj]). The symmetric square Sym 2 0( 1O) ... , ) = 
0(2,o, •••,()) © 0(o,o, •••,()) f° r so(2d) and d > 3 can be com- 
puted using Example 19.21 of Ref. (62|7 We verify the 
dimension of the commutant and show that (3) is a con- 
sequence of (2) by applying Proposition[9]which says that 
the dimension of the commutant of a representation is 
given by mf where the rrii are the multiplicities of the 
irreducible components of 0. For the rest of the proof we 
assume that condition (3) holds. We remark that the rep- 
resentation $|{ is irreducible as otherwise the dimension 
of the commutant would be larger than three. Thus, we 
obtain that t is semisimple. The dimension of the com- 
mutant allows only two possibilities: one of the restric- 
tions (Alt 2 <&)|t or (Sym 2 <I>)| t to the subalgebra t has to 
be irreducible. We emphasize that £ is given in an orthog- 
onal representation (i.e. a representation of real type) of 
even dimension, as t is given in an irreducible representa- 
tion obtained by restricting the standard representation 
of so (2c?). Therefore, we can use the list of all irreducible 
representations which are orthogonal or symplectic (i.e. 
of quaternionic type) and whose alternating or symmet- 
ric square is irreducible (Theorem 4.5 as well as Tables 7a 
and 7b of Ref. (6l|): (a) for su(2) the alternating square 
of the symplectic representation = (1) of dimension 
two, (b) for so(3) = su(2) the alternating square of the 
orthogonal representation = (2) of dimension three, (c) 
for so (2£ + 1) with I > 1 the alternating square of the 
orthogonal representation = (1, 0, • • • ,0) of dimension 
21+1, (d) for &o{2£) with £ > 3 the alternating square 
of the orthogonal representation = (1, 0, 0) of di- 
mension 2£, and (e) for sp(£) with £ > 1 the symmetric 
square of the symplectic representation = (1, 0, ■ ■ • , 0) 
of dimension 2£. Only possibility (d) fulfills all conditions 
which proves (1). □ 

Describing the matrices in the tensor square more ex- 
plicitly along the lines of Ref. [l~8| , we present a necessary 
and sufficient condition for full controllability in systems 
of quasifree fermions. 

Corollary 19. Consider a set of matrices {T„ | v G 
{0;l,...,m}} as given by Eq. (|22p generating the sys- 
tem algebra J C so (2d) with d > 3. We obtain t = so (2d) 
if and only if the joint commutant of {T„ <g> l 2 d + l 2 d <8> 
T v | v G {0; 1, . . • , m}} has dimension three. □ 

Along the lines of Eq. (|22|) , one can apply Corollarv[T9l 
to the matrices T corresponding to the generators of 
so (2d) of Theorem [15] For d > 3 one can verify that 
the commutant of the tensor square has dimension three. 
But for d = 2 one computes a dimension of four as 
so(4) = su(2) +su(2) is not simple. 



For illustration, note that two elements in the commu- 
tant are trivial, to wit the identity and the generator for 
the SWAP-operation between the two tensor copies. The 
third element does not yet occur in the unitary case de- 
scribed in [l8j]: it is the projector P$ onto the totally 
anti-symmetric state. To see this, recall that Ref. |63| 
implies that if the Hamiltonians {iH v \ v G {0; 1, . . . , m}} 
generate a system algebra of orthogonal type, then there 
is an operator S G SL(N) satisfying 

SHl + H V S = (24) 

jointly for all v G {0;l,...,m} as in [l8[. Using Kro- 
necker products and writing \S) := vec (S) [U, one sees 
that | S) is in the intersection of all the kernels of the 
tensor squares, so 

(JJ„g> 1 + 1® H V )\S) = |0) 
{H v <g> 1 + 1 ® H V )\S){S\ = 0^ 
<*\S)(S\{H V ®1 + 12>H U ) = N (25) 

and thus P s ■= \S) (S\ G (H v ® 1 + 1 <g> H v )' holds jointly 
for all v G {0; 1, . . . , m}; Oat denotes the zero matrix of 
degree N. 

VI. PURE-STATE CONTROLLABILITY FOR 
QUASIFREE SYSTEMS 

In this section, we present a straightforward criterion 
for pure-state controllability of quasifree fermionic sys- 
tems with d modes. We relate it to transitive actions on 
homogeneous spaces of type SU(2d)/U(d). 

A. Quasifree States 

A fermionic state p on J-'iC ) is called quasifree or 
Gaussian if it vanishes on odd monomials of Majorana 
operators and factorizes on even monomials into the Wick 
expansion form 

d 

tr(pm kl ...m k2d ) = ^sgn(Tr) tr(pm fejt(2p _ 1) Tnfe ir(2j)) ). 

7T p— 1 

Here the sum runs over all pairings of [1, ... , 2d] , i.e., over 
all permutations 7r of [1, . . . , 2d] satisfying 7r(2g — 1) < 
7r(2g) and Tr(2q— 1) < Tr(2q+ 1) for all q. The covariance 
matrix of p is defined as the 2d x 2d skew-symmetric 
matrix with real entries 

G% = i (Tr(pm p m q ) - 2,5p q ) . (26) 

Due to the Wick expansion property, a quasifree state is 
uniquely characterized by its covariance matrix. (Gen- 
eral references for this section include (l7l [65T - |68| .) The 
following proposition resumes a known result on these 
covariance matrices (see, e.g., Lemma 2.1 and Theorem 
2.3 in [66]), which will be useful in the later development: 
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Proposition 20. The singular values of the covariance 
matrix of a d—mode fermionic state must lie between 
and 1. Conversely, for any 2d x 2d skew- symmetric ma- 
trix G p with singular values between and 1 there exist 
a quasifree state that has G p as a covariance matrix. 



B. Orbits and Stabilizers of Quasifree States under 
the Action of Quadratic Hamiltonians 



The action of the time-evolution unitaries generated 
by quadratic Hamiltonians on quasifree states can be de- 
scribed by the following proposition (see Lemma 2.6 in 

My- 

Proposition 21. Consider a quadratic Hamiltonian 



H 



2d 
p,q=l 



i p m q ) 



which is defined by the skew- symmetric matrix T and acts 
on a quasifree state p a with a (skew- symmetric) covari- 
ance matrix G a . The time-evolved state (at unit time), 
p = e~ tH p a e tH is again a quasifree state with a (skew- 
symmetric) covariance matrix G b = OrG a 0^p, where 
O t :=er lT G SO(2d). 

Any skew-symmetric matrix G can be brought into its 
canonical form 



/On 
-v x 



OoGO'a 



\ 



-V2 







V 



v N 

-VN / 



using a (not necessarily unique) element Oq G SO(2g!) 
where {vi}f—\ denotes the singular values of G. This 
means that a quasifree state can be reached from another 
one by the action of quadratic Hamiltonians if their co- 
variance matrices share the same singular values (includ- 
ing multiplicities). Let us now recall another result re- 
lated to the singular values of the covariance matrices of 
pure quasifree states (Theorem 6.2 in [66[, and Lemma 1 
in [H): 

Proposition 22. A quasifree state p is pure iff the fol- 
lowing (equivalent) conditions hold for its covariance ma- 
trix G p : 

(a) The rows (and columns) of G p are real unit vectors 
which are pairwise orthogonal to each other. 

(b) The singular values of G p are all 1. 

Applying this result together with Proposition [5TJ we 
obtain the following theorem: 

Theorem 23. The set of quadratic Hamiltonians act 
transitively on pure quasifree states, and the correspond- 
ing stabilizer algebras are isomorphic to u(d). 



Proof. We have already shown that the singular values 
of the covariance matrices (with multiplicities) form a 
separating set of invariants for the orbits generated by 
quadratic Hamiltonians over the set of quasifree states. 
This means, according to Proposition that the pure 
quasifree states form a single orbit. 

As the set of quadratic Hamiltonians generate a transi- 
tive action over the pure quasifree states, the correspond- 
ing stabilizer subalgebras are isomorphic to each other. 
Consider a quadratic Hamiltonian H with the coefficient 
matrices A and B as given in Eq. (f2"Tj) and the Fock state 
which is the projection onto the Fock vacuum vector 
f2. pu is left invariant under the time evolution generated 
by H (pQ = e~ %Ht pne %Ht ) iff f2 is an eigenvector of H. 
We obtain that HQ = 



d i i 

Apqifpfq—fipq-^) + ~^Bpqfpfq ~ l^^pqfpfq 

',9=1 
d 



o 



P =i 



p<q 



By noting that and /J /Jfi (with p < q) are linearly 
independent vectors, we can conclude that a quadratic 
Hamiltonian H leaves the Fock vacuum invariant iff H = 

El,=i A Pq(flfqSpql)- In Theorem M of Sec. EE] we 
will show that these operators form a Lie algebra isomor- 
phic to u(d). □ 

Note that Theorem [23] identifies the space of pure 
quasifree states with the quotient space SO(2d)/U(d). 



C. Conditions for Quasifree Pure-State 
Controllability 

According to Proposition [231 a se t °f quasifree control 
Hamiltonians {Hi, . . . , H m } allows for quasifree pure- 
state controllability, if the corresponding Hamiltoni- 
ans generate the full quasifree system algebra, i.e. if 
(iHi, . . . , iH m )ue = so(2d). It is natural to ask whether 
this condition is also a necessary. Remarkably, it turns 
out that this is not case, which is shown by the following 
proposition: 

Lemma 24. Consider a quasifree fermionic system with 
d > 1 modes. Let K be the subgroup of SO (2c?) which is 
isomorphic to SO(2d— 1) and stabilizes the first coordi- 
nate; the corresponding Lie algebra is denoted by t. Then 

(a) via its adjoint action the group K acts transitively on 
the set of all skew- symmetric covariance matrices of pure 
quasifree states ( whose singular values are all 1 ); 

(b) the quasifree system is pure-state controllable if its 
system algebra is conjugate under SO(2rf) to {. 

Proof. We prove the statement (a) by showing that 
all pure quasifree states can be transformed under K- 
conjugation to the same pure state. We employ an induc- 
tion on d. The base case d = 2 can be directly verified. It 
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follows from Proposition l2"27b) that the skew-symmetric 
covariance matrix of a pure quasifree state can be written 
as 

<*-(-„*,)■ 

where v% denotes a normalized (2d— l)-dimensional vec- 
tor and A\ denotes a (2c?— 1) x (2c?— l)-dimensional skew- 
symmetric matrix. We consider the action of a general 
orthogonal transformation 1 © 0\ with 0\ G SO(d — 1): 

(i \( o v{\(i \_( o „{of \ 

^ Or) \- Vl Ajy 0\)~ {-O lVl O x A x O\) ' 

Since any (2d— l)-dimensional vector v\ with unit length 
can be transformed by an orthogonal transformation 
to (1, 0, 0, .., 0), we can choose 0\ such that v\0\ = 
(1,0,0, ..,0). We have (OiAiO\)n = as the trans- 
formed matrix is skew-symmetric. Again by Proposi- 
tion H^b) we obtain the transformed matrix as 

/0 1 \ 

-1 

o vly 

V -v 2 A 2 ) 

where v 2 is a 2d — 3 dimensional unit real vector and A 2 
is a (2c? — 3) x (2c? — 3) skew-symmetric matrix. Now the 
proof of (a) follows using the induction hypothesis. The 
statement (b) is a consequence of (a). □ 

We relate Lemma [M] to what is known about transi- 
tive actions on the coset space SO(2c?)/U(d). Only Lie 
groups isomorphic to SO(2c? — 1) and SO(2c?) can act 
transitively (i.e. in a pure-state controllable manner) on 
the homogeneous space SO(2c?)/U(c?) assuming d > 3. 
The case d > 4 is discussed in (69|. For d — 3 we have 
SO(6) S SU(4) and SU(4)/U(3) = CP 3 (where CP 3 de- 
notes the complex projective space in four dimensions), 
and it is known that only subgroups of SU(4) isomorphic 
to SU(4) or Sp[2) ^ SO(5) can act transitively on CP 3 
(see p. 168 of [7Cj or p. 68 of Q; refer also to (z3|). 

In most cases the so(k — l)-subalgebras of so(k) are 
conjugate to each other. More precisely, Lemma 7 of 
[Zl| states that for 3 < k {4, 8} all subalgebras of so(k) 
whose dimension is equal to (k— l)(fc— 2)/2 are conjugate 
to each other under the action of the group SO(fc). In 
particular, it follows in these cases that all subalgebras 
of so(k) with dimension (k — l)(k — 2)/2 are isomorphic 
to So(k — 1). Interestingly, the last statement holds also 
for k G {4,8} (see Lemma 3 of [z3); however not all of 
these subalgebras of so(fc) are conjugate. We obtain the 
following theorem providing a necessary and sufficient 
condition for full quasifree pure-state controllability in 
the case of d > 4 or d = 3 modes: 

Theorem 25. A quasifree fermionic system with d > 
4 or d = 3 modes is fully pure-state controllable iff its 
system algebra is isomorphic to SO (2d) or so(2d — 1). 



Proof. "=>": Proposition [231 identifies the space of pure 
states with the homogeneous space SO(2c?)/U(c?). As- 
suming d > 3 we summarized above that a group acting 
transitively on this homogeneous space is isomorphic ei- 
ther to SO(2c?) or SO(2c!- 1). Thus only the full quasifree 
system algebra so (2c?) or a system algebra isomorphic to 
so(2d — 1) can generate a transitive action on the space 
of pure quasifree states. 

"^=": As discussed, all so(2c?— 1) subalgebras are conju- 
gate to each other for d > 4 and d = 3. Lemma l2"4l b') 
then implies that any set of Hamiltonians generating a 
system algebra isomorphic to so (2d — 1) will allow for 
quasifree pure-state controllability. □ 

Note that the cases d = 2 and d = 4 are well-known 
pathological exceptions. The algebra so (4) breaks up into 
a direct sum of two so(3)-algebras which hence cannot be 
conjugate to each other. For d = 4, there are three classes 
of non-conjugate subalgebras of type so (7) in so (8) where 
two classes are given by irreducible embeddings and the 
third one is conjugate to the reducible standard embed- 
ding fixing the first coordinate (74| . 

On a more general level, Theorem [25] can be seen as 
a fermionic variant of the pure-state controllability crite- 
rion for spin systems [25| , which has been recently gener- 
alized from one-dimensional subspaces to fc-dimensional 
Grassmannians 1721 . 



VII. TRANSLATION-INVARIANT SYSTEMS 

Here we relate the respective system algebras for 
translation-invariant spin systems and fermionic systems 
with and without restriction to quadratic interactions as 
well as with and without reflection symmetry. Particu- 
lar attention is paid to which algebras are generated by 
Hamiltonians of restricted interaction length. Finally, a 
classification of the orbit structure of pure translation- 
invariant quasifree states is presented. This can readily 
be used to bound the scaling of the gap for a class of 
quadratic Hamiltonians that are translation-invariant. 



A. Translation-Invariant Spin Chains 

Consider a chain of L qubits with Hilbert space 
®iLiC 2 . The translation unitary Ut is defined by its 
action on the canonical basis vectors as 

U T \ni,n 2 , . . . , n L ) = \n L ,n 1: . . . (27) 

where rij G {0,1}. We will determine the translation- 
invariant system algebra which is defined as the maximal 
Lie algebra of skew-hermitian matrices commuting with 
the translation unitary Ut- 

Lemma 26. The translation unitary can be spectrally 
decomposed as Ut = ^2e = o exp(2Tri£/ L)Pe, and the rank 
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rg of the spectral projection Pg is given by the Fourier 
transform 



L-l 



re 



1 2^ L ^ exp(-2Trikl/L), (28) 



fe=0 



where gcd(L, k) denotes the greatest common divisor of 
L and k. 

Proof. The eigenvalues of Ut are limited to exp(2iri£ / L) 
with I £ {0, . . . , L— 1} as the order of Ut is L, i.e. U% = 
1. Hence, the corresponding spectral decomposition is 
given by Ut = Y^iZo exp(2Tri£/L)Pg. This induces a 
unitary representation Dt of the cyclic group %i, which 
maps the k-th power of the generator g £ of degree 
L to DT(g k ) = Uji- Note that Dt splits up into a direct 
sum Dt = ®<?e{o,...,L-i}(AO ffidim(P ' ) containing dim(P f ) 
copies of the one-dimensional representations satisfying 
D e {g k ) = exp(2mk£/L). Therefore, we determine the 
rank of a projection Pg by computing the multiplicity of 
Dg using the character scalar product 



L-l 



re 



= ^J2tv[D T (g k )} tv[De(g k )]* 

1 £_1 

= jYl tr i D T(g k )} exp(-2irik£/L). 



k=0 



The trace of DT(g k ) is equal to the number of basis vec- 
tors left invariant since DT(g k ) is a permutation matrix 
in the canonical basis. From elementary combinatorial 
theory we know that a bit string (m, n%, . . . ,til) is left 
invariant under a cyclic shift by k positions if and only if 
it is of the form 

(ni,n 2 , . . . , n gc d(i,fe), • . • , nj., n 2 , ■ ■ ■ , n gc d(i,fc)) • 

It follows that the number of [/^-invariant basis vec- 
tors and — hence — the trace of DT{g k ) — U^ is equal to 
2 gcd(L,fc)^ Thus ^ the multiplicities oiDg 

are given accord- 



ingly by rg = \ YXZl 2 gcd(L ' fc) exp(-2mk£/L). 



□ 



Note that a Hamiltonian commutes with Ut iff it com- 
mutes with all spectral projections Pg of Ut- Combining 
this fact with Theorem[10]we obtain a characterization of 
the system algebra for translation-invariant spin systems: 

Theorem 27. The set of translation-invariant Hamil- 
tonians acting on a L-qubit system generate the system 
algebra t(L) := s^^Tq 1 11(77)], where the numbers rg are 
defined in Eq. (|28l) . 



One can completely analogously show that for a chain 
consisting of L systems with N levels, the system alge- 
bra is equal to s[(Bg=Q u(rN,e)], where tn,£ denotes the 
Fourier transform of the function N gcd ( L > k \ 











® if L ~ M 


u T q 









B. Short-Ranged Spin-Chain Hamiltonians 

In most physical scenarios, we may only have direct 
control over translation-invariant Hamiltonians of limited 
interaction range M < L. The corresponding subalgebra 
tM(L), or \m for short (note that t(L) = ti(L)), is gen- 
erated by the operators 



L-l 
q=0 



for all combinations of Q p £ {l2,X, Y, Z} apart from the 
case when Q p = 1 2 for all p e {1, ... , M}. 

Building on [l6[ to a large part of this subsection, we 
analyze the properties of the algebras \m for different 
values of M, We generalize a key observation of [l6| 
(from the particular case when K divides L): 

Proposition 28. The trace of the product of U^ K with 
a tensor product of Pauli operators Qt £ {l2,X, Y, Z} 
can be computed as 









'L/c-1 


tr 


U T K (®QiJ 




Q(qK+p) mod L 






p=i 


q=0 



where c := gcd(if, L). 



(29) 



Proof. To simplify our calculations, let us introduce the 
notation v(l) = (K + £) mod L, note that (y o v)(£) = 
v(v(£)) = [2K + £) mod L, or more generally v°' p (£) = 
(pK + £) mod L. We can now write the action of U^ K 
on an arbitrary standard basis vector as 

U^ K |m, 722, . . . , ni, . . . , til) = 
\n v (i),n v (2), ■ ■ ■ ■ • , "•«(£))• 

Without loss of generality we can confine the discussion 
to the case where K < L. We complete the proof, by 
evaluating the trace in Eq. (|2"9"|) as 



tr 




tr 




Qi I Ut K 



(ni,---,n L \(<^}Qi\u^ K \n 1 ,...,n L ) 

n£{0,l} L \i=l ) 

^2 (n>u ■■■,n L \ I (£)Qi J |n„(i), . . .,n„(i)) 

ne{0,l} L \i=l / 



Yl( n e\Qe\ n v(e)) 

n€{0,l} L <=1 
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which can further simplified to 

L 

^(K+l) mod L/ 

ti6{0,1} l 1=1 

c L/c-1 

= E [ ( n v°"{p)\Qp\ n v°(<i+ 1 )(p)) 

n£{0,l} L P=l q=l 

~L/c-l 

Q[(qL/c+p) mod L] 

q=0 



ntr 

p=l 



□ 



This allows us to prove the following proposition: 

Proposition 29. Let M denote a divisor of L. Given 
the two elements HIm-i € t/\/-i oxid iHm G tut, we 
obtain that 

^{U^Hm-i) = and ti[(U qM - U T qM )H M ] = 0. 
hold for any positive integer q. 

Proof. The Lie algebras tm-i and i\j are generated by 
elements of the form of 



U T q and 




respectively. We consider all combinations of Q p £ 
{l2,X,Y, Z} apart from the case when Q p = I2 for all 
p e {1, ...,M-1} (or p e {1,...,M}). We introduce 
the notation F(a,W) := 



tr ( U a T qM i E U q T 

q>=0 




q'=0 



where a e {1,-1} and W £ {M, A/— 1}. Using Proposi- 
tion [28l we compute the formulas 

F(l,M-l) = iL Y[ ti[Q p ] = 0, 
P =i 

M-l 

F(1,M) = iLtr [QiQjif] J] tr[Q p ], 

p=2 
M-l 

F(-l, M) = iLtx [QmQx] ] [ tv[Q p }. 



It follows that the respective statements in the proposi- 
tion hold for the generators of Im-i and Im- Now we 
prove this consequence also for any element in Im— 1 (or 
tj\/). First, let us note that the elements generated must 
be contained in [Im-i, tM-i] (or [Im^m])- Second, since 
all elements in commute with [Zf> , we have that 
tr(U!^ M ih) = holds for any element ih € [tjv/,t/\/], as 

tT{[ih l M ,ih M ]Uf) 

= tr(i/i^ ih 2 M f/j. M ) — ti(ih M ih M J7y ) 



= tr(ih\j ih M U!f M ) — ti(ih M U!f M ih M ) 

= tT(ih M ih 2 M U* M ) - tr{ih l M ih 2 M U^') = 0. 



□ 



Theorem 30. Let AI denote a proper divisor of the num- 
ber of spins L, and define R := L/M. We obtain: 

(a) The algebra t/\/-i is isomorphic to a Lie subalgebra 

of [J2e=o su ( r ^ )] + Ei=i u (l)] an d does not generate ij, . 

(b) The algebra tjvj is isomorphic to a Lie subalgebra of 

[J2e=o su ( r t)] + Ei^i + u (l)] an d does not gener- 
ate {l- 

(c) iM-xi 1 hi- 

Proof, (a) Obviously for any integer q the following holds 

L-l L-l 



U, 



qM _ 



E exp(2mqM£/L)Pi = ^ exp(2niq£/R)P e 



R-l 



t=0 



jqM 



= E exp(2mq£'/R) ( £ P pR+ t J . 

l'=0 \p=0 / 

One can invert these equations as 

(M-l \ R-l 

E P PR+C =nE exp(-27riq£'/R) U q 
p=o ) 9=0 

If ih G t.M_i, we obtain by applying Proposition 1291 that 

tr (ih p pR+A = 0- (30) 
V p=q J 

The maximal Lie algebra consisting of skew-hermitian 
matrices which satisfy the condition in Eq. ()30[) is iso- 
morphic to ®^zl(s[® p ^r 1 u(r pR+il )}) - Efco 1 SU M] + 

EiLi^^l)]- Hence tjw-i must be a subalgebra. 

(b) For elements ig <G iu , Proposition [29] and Eq. (|30|) 

imply that 



tr 



i.9 E ( P pR+e' ~ PpB+L-t) 
p=o 



= 0. 



(31) 



The maximal Lie algebra consisting of skew-hermitian 
matrices which satisfy the condition in Eq. (pH]) is iso- 



morphic to [£^ =0 su(r £ )] + E i= i 
(c) Let 



u(l)]. 



L-l 



p=2 



i/i = iJ2 U T [ x ® if M " 2 ® ^ ® 1 2 

<j'=0 
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Obviously, ih G tjvf. Using the formula for F(+l, M— 1), 
we obtain tr(z/i) = iL2 L ~ 1 . Hence, ih ^ Im-i- D 

In particular, this theorem implies that the algebra tt 
of translation-invariant Hamiltonians cannot be gener- 
ated from the subclass of nearest-neighbor Hamiltonians, 
cp. also [l6[. More precisely, one finds: 

Corollary 31. // L is even, ti is isomorphic to a Lie 
subalgebra of [J2e=o su ( r ^)] + E*=i u (l)]- F° r °dd L, 
ii is isomorphic to a Lie subalgebra of E_j_ su(r^)] + 



E 



(L-l)/2 



u(l)]. 



rjr of the spectral projection Pp is given by the Fourier 
transform 

^ d-l 

rf := - h± (d, k) exp(-2Ttik£/d) (34) 

fc=0 

ofh±(d, k) where I G {0, . . . , d — 1}, 



h+(d, k) 



2 gcd(d,fe) y y gcd ( d) fc ) 

is even, 

2 gcd(d,fe)-i ifd/gcd{d,k) is odd, 



and 



C. Translation-Invariant Fermionic Systems 

To determine the system algebra generated by all 
translation-invariant Hamiltonians of a fermionic chain, 
we can follow similar lines as in Sec. IVII A[ Here, how- 
ever, we additionally have to consider the parity su- 
perselection rule. We define the fermionic translation- 
invariant system algebra as the maximal Lie subalge- 
bra of su(2 d_1 ) ©su(2 d_1 ) [see Theorem Q] which con- 
tains only skew-hermitian matrices commuting with the 
fermionic translation unitary IA, which is defined below 
such that it commutes with the parity operator P (see 
Eq. (H3|). The standard orthonormal basis in the Fock 
space for a chain of d fermionic modes is given by 

\ ni ,n 2 ,...n d ) :=(f!) n Htir |0> (32) 

with Hi G {0, 1}. Note that for the purpose of unambigu- 
ously defining this basis, we order the operators (f}) ni in 
Eq. (|32f with respect to their site index i. The fermionic 
translation unitary Li is defined by its action 

U |m,n 2) ...n d )=U (fir UlT 2 ■ ■ ■ (fl) nd |0) 

= (/ 2 t r---(/i) nd - i (/i t ) rid io> 

= (— i)™ d (™ i +" 2 +---+ n < i - 1 ) (fl) ,ld ( f 2 ) ni • • • (f d ) nd ~ 1 |o) 
= (_l)»'(»i+»»+-+'w-0 \ nd , ni ,..., n d _i) (33) 

on the standard basis. The adjoint action of U on the 
creation operators f\ is then given by 

U f\lfi = / ( ^ +lmod d y 

The superselection rule for fermions splits the spectral 
decomposition of the translation unitary into two blocks 
corresponding to the positive and negative parity sub- 
space. The translation unitary U commutes with the 
parity operator P, and hence U = U+ + U- is block- 
diagonal in the eigenbasis of P where U+ := P+UP+ and 
U- := P-UP-. The following Lemma gives the spectral 
decomposition of the operators IA± : 

Lemma 32. The unitary operators IA± can be spectrally 
decomposed as IA± = E^=o e 27!li / d P^ , where the rank 



h-(d, k) 




if dj gcd(d, k) is even, 
if dj gcd(c?, k) is odd. 



Proof. We determine the spectral decomposition of U+ 
and W_ along the lines of Lemma [2^1 Let T + (C d ) denote 
the subspace spanned by those basis vectors of Eq. (f3"2"f 
for which n = Ei=i n i 1S even - Likewise, J r _(C d ) cor- 



responds to the case of odd n. 



As Ui 



1 



the eigenvalues of IA± are of the form exp(27rz£/d) with 
i G {0, . . . , d— 1}. Hence, the spectral decomposition is 
given by U± = Y^eZo Gxp(2iri£/d)Pp. We define rep- 
resentations T>± of the cyclic group 7j d which map the 
fc-th power of the generator g G 1 d of degree d to 
l- ) ±(9 k ) '■= U±. Note that T>± splits up into a direct sum 
V± = ®ee{o,...,L-i}{De) (Bdiln{Pe} containing dim(P/) 
copies of the one-dimensional representations satisfying 
Di(g k ) = exp(2irik£ / d) . The rank of the projection 
Pp is equal to the multiplicity of Di in the decomposition 
of the reducible representation T>± . This multiplicity can 
be computed as the character scalar product 

d-l 



Wtr[D ± (g k )} tr{D e (g k )}* 



k=0 
d-l 



^5>[2? ± ( 5 fc )] exp(-2mlk/<t). 



k=0 

In the standard basis, all matrix entries of T>±(g k ) 
are contained in the set {0,1,-1}. Moreover, exactly 
one entry in each row and column differs from zero. 
The diagonal contains — as a consequence of Eq. ([33)) 
only elements from the set {0,1}. Hence the trace of 
T> + (g k ) and of T>-(g k ) equals the number of W fe -invariant 
\ni,n 2 , ■ ■ .n d ) basis vectors with the n = Ej=i n i 
ing even and odd, respectively. Recall from the proof of 
Lemma 1261 that a bit string (m, n 2 , ■ ■ ■ , tin) is left in- 
variant under a cyclic shift by k positions iff it is of the 
form 

(ni,n 2 , . . . ,n gc d(Ufc), ■ ■ ■ ,ni,n 2 , . . . , n gc d(<ifc)) • 

If dj gcd(<i, k) is even, we obtain that the sum n is even 
for all of the 2 gcd ( d,fc ) bit strings invariant under a cylic 
shift by k positions; thus it follows that the traces of U\ 
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and are 2 gcd ( d,fc ) and zero, respectively. If dj gcd(d, k) 
is odd, the sum n is odd for one half of the 2 gcd ( d,fc ) bit 
strings and even for the other one; hence both traces 
tr(W±) are equal to 2 gcd ( d,fc ) _1 . This completes the proof. 

□ 

The above Lemma together with Theorem [TU] implies 
the following characterization of the system algebra for 
translation-invariant fermionic systems: 

Theorem 33. The set of translation-invariant Hamilto- 
nians of a d-mode fermionic system generate the system 
algebra s[®^~q u(r^)](Bs[®g~Q u(r^ )], where the numbers 



£l,£f, andij are lin- 



are defined in Eq. ([34] 



D. Translation-Invariant Quadratic Hamiltonians 

A quadratic Hamiltonian H is translation-invariant 
(i.e. [H, U] = 0) iff the coefficient matrices A and B 
in Eq. (HU are cyclic (i.e. A nm - A n+ i^ m+1 = B nm - 
B n +i,m+i = 0). To study such Hamiltonians, it is use- 
ful to rewrite them in terms of the Fourier-transformed 
annihilation and creation operators 

f k :=±f: f P e-^ d and ft := -L £ 



P =i 



Vd 



p=l 



with k e {0, 1, . . . , d— 1}, which satisfy again the canoni- 
cal anticommutation relations 



{flfl}={f k J k '} = ^d {flf k ,} = 5 kk ,l. 



(35) 



A Hamiltonian from Eq. (|2lj) with cyclic A and B can 
now be rewritten as 

d-1 

H = E Mflf k !) + \B k flf\_ k (36) 

k=Q 

applying the definitions A^ :— Yl P =i Ai p exp(— 2iripk/d) 
and Bk '■= 2~2 p= i Bi p exp(—2iripk/d), as well as the no- 
tation fa = fo. The hermicity of A and the skew- 
symmetry of B translates into the properties A k = A* d _ k 

and Bk = —Bd-k- This allows us to decompose the 
Hamiltonian into a four-part sum 



L(d-1)/2J L(d-i)/2j 

k=l k=l 
L(d-1)/2J \_d/2\ 



fe=i 



fc=0 



where t\ := i{flf k -fl_J d _ k ), % ■= (fUL k +f d -J k ), 
?J ■■= i(flfl k -f d -J k ), aud^| := (fVk+flkf d -k-^) 
with k G {1, . . . , L(d- 1)/2J} as well as £^ /2 = (/] /2 / d/2 - 
1/2) (for d even) and £^ = (f£f - 1/2). Note that the 



operators £^ 2 (for d even), £%, ££, 
early independent and span the ( [d — lj + d)-dimensional 
space of all translation-invariant quadratic Hamiltonians. 
For notational convenience we also introduce the dummy 
operators £^ 2 = (assuming d is even) and £^ = for 
Q G {1,X, Y}. We can now characterize the system al- 
gebra: 

Theorem 34. The system algebra q d of all translation- 
invariant quadratic Hamiltonians with d modes is iso- 
morphic to the Lie algebra u(2)] +u(l) for odd 
d and to [J2[t~ 1 2)/2 u(2)] + u(l) + u(l) for even d. 

Proof. If d = 2m — 1 is odd, the generators i£\ 



k- 



k ■ 



ilY , i£|, and Hq can be partitioned into m pairwise- 
commuting sets which each generate linear subspaces as 

L :={il%) m and L k := {i£l,Uf,H k 

with k G {1, . . . , m— 1}. The commutation proper- 
ties [L k ,L k >] = (with k ^ k') follow from Eq. (|55)) . 
Moreover, Lq is one-dimensional and forms a u(l)- 
algebra. Using Eq. (|35|) . the relations [flf a , flfy] = 

([/l/aJa/J) 1 = flfl and [/t/7,/ 6 / ] = /l/ a + / 6 t /6-l 
can be deduced for a 7^ b. Substituting k and d — k into 
a and 6 in the previous formula, one can verify directly 
that the correspondence 

ill h- it, ilf 1 — ^ iX, ^ i — ^ iY, ^ i-> iZ 

provides an explicit Lie isomorphism between L k and 
u(2). If d = 2m is even, the system algebra consists of 
the above-described generators supplemented with the el- 
ement i£^ 2 - This additional element commutes with all 
the other generators and — therefore — provides an addi- 
tional u(l). □ 

The isomorphism between L k and u(2) as given in 
the proof leads to a compact formula for the time evo- 
lution (in the Heisenberg picture) of the elements of 
Lk- Since the operators , £\, £f, and £\ (with k G 
{1, . . . , [(d — 1)/2J }) satisfy the same commutation rela- 
tions as the Pauli matrices X, Y, Z, and 1, their time- 
evolution generated by the Hamiltonian H in Eq. (|37[) can 
be straightforwardly related to a qubit time-evolution 



JHt- 



-H (ai4 + ax£f + a Y £j + az£ k ) 



z^ e -im 



ia-ill 



JG{X,Y,Z} 



ia Q £% ti(e iH 'Qe- iH -Q), 



(38) 



where H s = <Ke(A fe )Z + 5Ke(E fe )X/2 + Jm(B k )Y /2. 

The twisted reflection symmetry plays an important 
role in translation-invariant quasifree fermionic systems. 
It is defined by the unitary 

n\ nil n 2 ,...n d ) =i^^ n ^\n d ,n d ^ ll ...n 1 ). (39) 
whose adjoint action on creation operators is specified as 



(d-£+lmod d)' 
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A given translation-invariant quasifree Hamiltonian is 1Z- 
symmetric (i.e. [1Z, H] = 0) iff the coefficient matrix is 
restricted to be real. In our language, these Hamiltonians 
are exactly the ones for which 3m(Ak) = 0, i.e., they 
are spanned by the operators i£^/ 2 (f° r ^ even), Hq, il^, 
ii^, and ii^. From the proof of Theorem [34] one can 
immediately deduce the corresponding system algebra: 



Corollary 35. Let denote the minimal system al- 
gebra which contains all d—mode fermionic Hamilto- 
nians that are quasifree, translation-invariant, and 1Z- 
symmetric. The system algebra is isomorphic to 
E^ 1)/2 su(2)]+u(l) /or odddand to [Y}^' 2 Sll(2)] + 
u(l) + u(l) for even d. 



Given the system algebras q^ and qj, we investigate 
the subalgebras generated by short-range Hamiltonians. 
It will be useful to introduce for p G {1, . . . , [(d ~ 1)/2J } 
the Hamiltonians 



metric identities, we obtain 



L(d-1)/2J 



(40a) 



k=0 



L(d-l)/2j 



X 

k i 



(40b) 



k=0 



L(d-l)/2j 



Y 
A- • 



e=i 



k=0 



(40c) 



Ld/2J 



k=0 



(40d) 



as well as the additional ones (h^, 2 only for d even) 



Ld/2J 



ho ■= \ Y^Ulft +f\h - 1 ) = E £ t (40e) 

i=l k=Q 

h z d/2 ■= \ Y.{flh +p +flvh ) - E(-i) ft 4 (4Qf) 



[ihl,ih%] =[ihl,ih£] = [ihl,ihZ]=o, 



\ih*,ihj] = 



\ihl,ihf) 



Kfya+ft) mod |_d/2J mod |_d/2J )' 



|(sgn(d-a-6) hf a+b) mod [d/2j 
- sgn(a-6) hf a _ b) mod Ld/2J ) 

|(sgn(d-a-&) hj a+b) mod Ld/2J 
-sgn(a-6) hJ a _ b)mod [d/2j ) 



(41a) 

j). 
(41b) 

(41c) 
(41d) 



for a, b G {0, . . . , Ld/2J}. In [16[ it was shown that al- 
ready the nearest-neighbor Hamiltonians of q J generate 
the whole qj. We are now in the position to provide a 
more systematic proof of their result: 

Proposition 36. The system algebra qj can be gen- 
erated using an one-site-local Hamiltonian ih$ and a 
nearest-neighbor Hamiltonian i{a\hf-\-a 2 h^+a 3 hf) with 
a.i G M assuming that a 2 ^ or a?, 7^ /or odd d and 
additionally requiring ai 7^ /or euen 0!. 

Proo/. (1) From Eq. (gl]) we know that i/ig, ihf, ihf, 
and i/i^ would generate the whole q^. (2) Suppose that 
ai 7^ and a 2 + a\ ^ 0. From 2[i/ig , i(a\hf + a 2 hf + 
a 3hJ)] = cn 2 ihj — a^ihf and 2[ihQ,a 2 ihJ — a^ihf] = 
—a 2 ihf — a^ihj it follows that one can generate ih$, 
ihf, ihf, and ihj . Hence according to observation (1), 
the whole qf is generated. (3) Suppose now that ct\ = 0, 
d is odd, and a 2 + al ^ 0. From 2[ih^ , i(a 2 tif + a 3 hj)] = 
a 2 ihj — a^ihf one can generate ih$ , ihf , and ihf . From 
Eq. (j4"Tj) it follows that these generators in turn generate 
all ih 2pmod d . Since d is odd, ihf is also generated. Hence 
we obtain ih%, ihf, ihf, and ihf , and according to (1), 
the algebra q J is generated. □ 

For the more general q^, we obtain a slightly larger 
system algebra when we do not assume 7?.-symmetry: 

Lemma 37. The Hamiltonians in q^ with interaction 
length smaller than M (with M < \_d/2\) generate a 
system algebra which is isomorphic to Eiii su (2)] + 
EtLi u(l) for odd d and to \£<£»' 3 su(2)] + u(l) 
for even d. 

Proof. From Proposition [55] we know that the operators 
with Q G {X, Y, Z} already generate q J which is 

isomorphic to Ei=i ^(2)] + u (l) f° r °dd d and to 
E,=T 2)/2 su(2)] +u(l) +u(l) for even d. We have M-l 
additional operators /jj with q G {1, . . . , M— 1} which are 
linearly independent and commuting. These generate the 
other parts corresponding Ef=i 1 u (l)- I— ' 



In these definition we used cyclic indices, e.g. fd+a = fa- 
The operators h^, 2 (for d even), h%, hp, hp, h£, and hj 
span q^ linearly. Using the identities above, the commu- 
tation relations of the £j} operators, and some trigono- 



E. Orbits of Pure Translation-Invariant Quasifree 
States 

We characterize now the orbits of pure translation- 
invariant quasifree states under the action of translation- 
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invariant quadratic Hamiltonians. Since the opera- 
tors l\ = i(flh~fd-kfd-k) commute with all the 
other translation-invariant quadratic Hamiltonians (as 
discussed in Sec. IVIIDp . their expectation values stay 
invariant under the considered time evolutions. At the 
end of the section, we show that these invariant expec- 
tation values even form a separating set of invariants for 
the orbits of pure translation-invariant quasifree states. 

Let us recall that a quasifree state is fully characterized 
by its Majorana covariance matrix, defined in Eq. (|26|) . 
The translation unitary U acts on the Majorana oper- 
ators by conjugation as U m p W = m( p+2mo( j 2d)- It 
follows that a quasifree state p is translation-invariant 
(i.e. [p,U] = 0) iff its covariance matrix G vq is doubly- 
cyclic, i.e. G P q = G( p+2mo d 2d),( 9 +2mod 2d)- The double- 
cyclicity of G implies that it can be expressed as a block- 
Fourier transform of a block-diagonal matrix, i.e. 



with any translation-invariant Hamiltonian, we have that 



G = Uf GUp , 



(42) 



where U F := (J ?) ® W with W pq := exp(27ri/d) q - p and 
G = ©felj«g(fc) with g(k) being 2 x 2-matrices. The ma- 
trices g(k) can be calculated by the inverse block-Fourier 
transform 



2-rrkti/d 



G\ : 2l-1 Gi t 2t 
Gl,U-X G2,2£ 



(43) 



The fact that G is skew-symmetric and real implies 

g(d-k) = ~~g T (k). (44) 



Moreover, due to Eq. (|4"2"|) the set of eigenvalues of all the 
matrices g(k) equals the one of —iG (including multiplic- 
ities). Combining these observations with Proposition |2"01 
and Proposition [22] we obtain the following characteri- 
zation of pure translation-invariant quasifree states: 

Lemma 38. A set of 2 x 2 matrices g{k) (with k £ 
{0, . . . , d— 1}) defines a covariance matrix of a pure 
quasifree state through Eq. (|42p iff they satisfy Eq. (|44p 
and their eigenvalues are in the set {1,-1}. 

The entries of g(k) and the expectation values of the 
Lfc operators can be related by 



is(fc) = l 2 (4)+X<#)+Y<#)+Z<^) 



(45) 



using Eq. 03J and the definitions for t\, if, £*, and £%. 
This allows for proving the main theorem of this Section: 

Theorem 39. Two pure quasifree states p\ and p2 can 
be connected through the action of a translation-invariant 
quadratic Hamiltonian if and only i/tr(pifl) = tr(p2^J) 
holds for all l\ with k £ {0, . . . , |_(d-l)/2j}. 

Proof. First, we consider the 'if-case: Let H be a 
translation-invariant quadratic Hamiltonians for which 
pi = e~ lHt p2e lHt holds. Since the operators l\ commute 



tr(e- iHt p 2 e im £l) = 



tr(p 2 e lHt e 



Htol„-iHt 



tr(p2^fe)- Second, we treat the 'only if'-case: Let gi(k) 
and 02 (k) denote the Fourier-transformed Majorana two- 
point functions (defined as in Eq. (|43[) ) of p\ and p 2l re- 
spectively. The action of a translation-invariant Hamil- 
tonian, p a i y e~ lH p a e lH is represented by the map 



g a (k)^U(k)g a (k)U(k^ 



(46) 



where U(k) is given by exp[-i9le(A fe )Z - mt(B k )X/2 - 
Om(JBfc)Y/2]. Using Eq. 05]), we obtain tr(p Q ^) = 
itr[g a (k)] for a £ {1,2}. These expectation values have 
to be in the set {—2, 0, 2}, since the eigenvalues of g\{k) 
and g2(k) are in the set {—1, 1}. Then, it follows from 
tv{pxi\) = tv{p2(-}.) that the expectation values of gi(k) 
and §2 (k) coincide. Thus, we obtain from Eq. (|4T>|) that 
pi and p2 can be transformed into each other. □ 

Finally, we turn to the 7?.-symmetric setting, as intro- 
duced in Sec. I VII Dl and determine the orbit structure of 
quasifree pure states which are translation-invariant and 
Asymmetric under the action of Hamiltonians in qj. 

Proposition 40. The Hamiltonians from the algebra q J 
act transitively on the set of quasifree pure states which 
are translation-invariant and IZ-symmetric. 

Proof. Since Iti^Rr 1 = the expectation value of 

these operators in 7?.-symmetric states must vanish as 
tr(p^) = -tT{p1Upl- 1 ) = -tx{n- l P 1Ul) = -tx{pll). 
Moreover, by Theorem [35] we know that two pure 
translation-invariant states are on the same q^-orbit iff 
the expectation values of the i\ operators coincide for all 
k £ {0, . . . , \_{d — 1)/2J}. Hence the translation-invariant 
1Z— symmetric states lie on the same q^-orbit. As Eq. (|4*6"]1 
implies that the q^-orbits are equivalent to q J-orbits, wc 
have proved the Proposition. □ 



F. An Application to Many-Body Physics 

In many-body physics, one of the important character- 
istics of quantum criticality is the closing of the gap. This 
means that the energy difference between the ground 
state and the first excited state goes to zero in the ther- 
modynamic limit, when the number of spins or fermionic 
modes goes to infinity. Quasifree fermionic models can 
display both gapped and gapless behavior. Using the 
techniques developed in the previous subsections, we 
will prove that the gap always disappears (closes) for 
translation-invariant quasifree models if the coefficient 
matrix A of Eq. (|21| is purely imaginary while B is an 
arbitrary, complex skew-symmetric matrix. 

To formalize this statement, let us consider a set a r 
of fixed (finite) real numbers with r £ {0, . . . , M} and a 
set b r of fixed complex numbers (of finite modulus) with 
r £ {0, ...,M}. With these stipulations, we define for 
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any d > 2M the cylic d x d matrices Ad and Bd (or A 
and B for short) by specifying their entries 



another quasifree state p d which is defined through its 
Majorana two-point functions 



4 — 

^pq ■ — 



-ia v - 



\fq- P e{0, 
if P -qe{0, 

otherwise, 



,M}, 
,M}, 



(47) 



and 



B 



pg 



Jq-p 



iiq-pe{l,...,M}, 
i£p-q€{l,...,M}, 
otherwise. 



(48) 



By applying these definitions to Eq. ([21]) we obtain: 

Theorem 41. Given the positive integers d and M 
with d > 2M , consider the corresponding translation- 
invariant quasifree Hamiltonian 

d 

Hd = / ] Apq^f^fq — Spq^+^Bpqf^f^—^Bpqfpfq, 
p,q=l 



where A andB are defined in Eqs. (|47| and (|48|l . Assume 
that Hd has a unique ground state. The gap Ad of Hd 

is bounded as Ad < 8?r rf M X)p=i(l a pl + I^pDj *- e - ^ e ff a P 
closes algebraically in the thermodynamic limit of d going 
to infinity. 



Proof. Since Hd is translation-invariant and its coef- 
ficient matrix is imaginary, it can be decomposed in 
terms of the operators with Q e {1,X,Y} and 
fcG{l,...,L(d-l)/2j}afl 



L(d-1)/2J 

H d = J2 ~ a ^l 

fe=i 



\Pk^k + h^k^k ) 



using the definitions := — X^pli a p 
5f == -«e[E^i& P sin(-27rpfc/d)] 
— ^ m Epli sin(— 27rpfc/eJ)]. Let pd be a pure quasifree 
state, and let ga(k) denote its Fourier-transformed Ma- 
jorana two-point functions (see Eq. 143]) . From Eq. 
we know that pd is an eigenstate of Hd iff [fr/f X + 
b^Y,gd(k)] = 0. The eigenvalue of Hd corresponding 
to this state is given by 



-2irpk/d), 



as well as b\ := 



L(d-1)/2J 



tT:{pdH d ) = 2J tr[ig d (fc)(a fe l 2 - 



^X+I^Y)]. 



fc=i 



(49) 



Let us emphasize that the proof builds on the fact that 
M is fixed and finite, while d goes to infinity in the ther- 
modynamic limit. Among the eigenstates of Hd, consider 
the (unique) ground state p d s , whose Fourier-transformed 
Majorana two-point functions (see Eq. |4"5|) will be de- 
noted by g d s (k). From this ground state let us construct 



ffe(l) 



12 if 9 d gs(X)± 

-It otherwise, 



while for general k ^ 1 we assign g d (k) := g d s (k). 

The pure quasifree state p d is an eigenstate of Hd, 
since according to Eq. (|46)) its Fourier-transformed Majo- 
rana two-point function stays invariant during the time- 
evolution generated by Hd- Using Eq. (|49| . we can cal- 
culate the difference between the energies corresponding 
to Pg S and p d as 

A d := ti[(p d e - p d gs )Hd] 

L(d-l)/2j 

= E tr ([ffeW - ~9 d gs (k)}(a k l2+fb§X+fblY) 
k=l 



tr 



([3f(l)-5 3 d s (l)](5il 2 +i6fX+i6rY) 



<2||[gf(fc)-g 9 rf s (A : )](a fc l 2 +pfX+prY)|| 



\bi 



<4(|5i 

M 

a p sin(27rp/<i) 



< 4 



p=i 



6 P sin(27rp/ ci) 



M 



< — 2_> p i + i&pI) 



This completes the proof of the theorem. 



□ 



VIII. PARTICLE-NUMBER CONSERVING 
SYSTEMS 

Finally, we treat fermonic systems with fixed particle 
number. Their system algebras are given both for the 
general case as well as for the quasifree case. 



A. The System Algebra of Particle-Number 
Conserving Hamiltonians 

Let P n denote the orthogonal projection from the Fock 
space J"(C d ) = ©f l=0 A™ C d onto the n-particle sub- 
space A n C d C T(C d ) of dimension ( d ). The particle- 
number operator n of a fermionic system is defined as 
n := J2 n =o nP n- Note that Y, p =i fpfp^n = nip n holds 
for any ip n € A n C d . Hence, the particle number operator 
can also be expressed as 



fpfp 

p=i 
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A fermionic Hamiltonian H is called particle-number con- 
serving if it commutes with n. Using Theorem 1101 one 
directly obtains the corresponding system algebra. 

Proposition 42. The system algebra of particle-number 
conserving fermionic system with d modes is given by 



TABLE I. System algebras for d-mode fermionic systems 




e u[cd 



: odd 



B. Quadratic Hamiltonians 



A quadratic Hamiltonian H is particle-number con- 
serving iff its coefficient matrix B of Eq. (f2"Tj) is zero, i.e., 

iff H = Y. d p , q =i A pq(flfq ~ <Wf) where A denotes any 
hermitian matrix. The corresponding system algebra is 
given by following Proposition. 



Proposition 43. The system algebra of the particle- 
number conserving quadratic d-mode Hamiltonians is iso- 
morphic to u(d). 

Proof. Let l denote the R-linear mapping from the d- 
mode Hamiltonians which are quadratic and particle- 
number conserving to the d x d skew-hermitian matri- 
ces. We define l using t(i(flf q +flf p -^)) = i(e pq +e qp ) 



and i(f I f 9 -flf p ) 
with entries [e pq ] uv 



-e qp , where e pq denotes a matrix 
$ P ud q v Note that the canconical 



anticommutation relations imply that [f p f q ,frf s 

8 pS f q f}+8qrfpf s -SpsSqA-5psS qr 8 pq (f q f}.+f P ~f s -^)- 

Thus, l is a homomorphism since 

i{ [ K± {flf q ±f\f p -S pq |), k ± (/t/ s ±/t/ s -5 rs f )] ) 
= K 2 ± i([5 qr (flf s Tf f s f p ) ± S pr (flf s Tflf q ), 

5 P s(flf r Tftf q ) ± SqsifUrTfUp)]) 

= \k± (e pq ie qp ) , k± (e rs ie sr )] 

= [t(«± {flf q ±f}f p -6 pq ±)), l(k ± {ftf s ±ftf s -6 rs I))], 

where n + = i and K- = 1. The map l is even an isomor- 
phism as its kernel is trivial. The proposition follows as 
the Lie algebra u(d) is isomorphic to the Lie algebra of 
d x d skew-hermitian matrices. □ 



Symmetries' 1 



System algebra Details 



general systems: su(2 d_1 ) su(2 d_1 ) Thm | 

{T} s[® d e Zo u(r+)] ©s[®£oU(r7)] Thm. [33] 

{N} s(e„cvcnu[Q])ffis(® nodd u[( d )]) Prop.lH 



quasifree systems: 
{T},d odd 
{T}, d even 
{T,R},d odd 
{T,R},d even 
{N} 



so(2d)^ Prop. [13] 
[E^T 1,/2 w(2)]+u(l) Thm. EH 
[£^7 2)/2 "(2)]+u(l)+u(l) Thm. [34] 
E£ 1)/2 su(2)]+u(l) Cor. [35] 
[££ 2>/2 *u(2)]+u(l)+u(l) Cor. [35] 
u(d) Prop. [43] 



a besides parity superselection rule P: T = translation-invariance, 
R = twisted reflection symmetry, N = particle-no. conservation 

b the orthogonal algebra is represented as direct sum of two equal 
copies given as irreducible blocks of dimension 2 , the system 
algebra so(2d) itself was determined already, e.g., in Ref. l3oll . 



interactions (quasifree particles) as well as with and with- 
out symmetries beyond parity such as translation in- 
variance, twisted reflection symmetry, or particle-number 
conservation. Once having established the system alge- 
bras, the group orbits of a given (pure or mixed) initial 
quantum state determine the respective reachable sets of 
all states a system can be driven into by coherent con- 
trol. The special case of pure-state reachability and its 
relation to symmetric spaces has been treated with par- 
ticular attention. On a general scale, the system algebras 
determined serve as a dynamic fingerprint. 

Their application to quantum simulation has been elu- 
cidated in a plethora of paradigmatic correspondences 
between fermionic systems and spin systems. Hence we 
anticipate the comprehensive findings presented here will 
find a broad scope of use. 
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Appendix A: Parametrizations of Quadratic 
Hamiltonians 



In this appendix, we discuss various parametrizations 
of quadratic Hamiltonians related to the one of Eq. (|2"Tj) 
in Sec. [Vj We start with the parametrization 



with W„ 



H 



p.q 



=1 



fpfq + Dpqfpfq + E pq f p f q + F pq f p f q 



by complex d x d-matrices C , D, E, and F. Hermiticity 
of H requires C = F+, D = D\ and E = E\ while the 
(anti-) commutator relations enforce C = — C*, D = — £7*, 
and F = —F*. Setting A := 2E and B := -2C*, we 
recover the notation of Eq. (j2"Tj) and obtain 



1 d 

H = 2 ^ ] ~Bp q fpfq — Ap q fpfq + A pq fpf q + B p qf p f q 
p,q=l 
d 



*Rt(B pq )(l- IXZ- ZXI- I-I-IYZ- ZYI- I) 

g— 1 p— g— 1 d— p q—l p— g— 1 d— p 
q—l p—q—1 d—p q—l p—q—1 d—p 

3m(B pg )(l^YZ^XT^+l^XZ^YI^a) 

q—l p—q—1 d—p q—l p—q—1 d—p 

Jm(A P9 )(l--IYZ- - ZXI - - I-I--IXZ- - ZYI - - I). 

q—l p—q—1 d—p q—l p—q—1 d—p 



Appendix B: Applications of Quasifree Fermions to 
Spin Systems 

Here we take new fermionic approaches to exhaustively 
prove and improve some results of Ref. [l8| , where some 
proofs were still sketchy — thereby also filling a desidera- 
tum voiced in 1751. 



1 T3* f f i o A ( ft f x i\ i d ft ft I - A Spin System with System Algebra $o(2n + 1) 

9 Pq JPJq ' 2J\ p q \ Jp J q ®pq 2 ' Pq J p Jq 



p,q=i 



\ £ Kt{B pq )(flfl-f p f q )+Kt{A m )(flf q -f p f t 



Proposition 44 (see Proposition 27 in [18|])- Consider 
a Heisenberg-XX. chain with the drift Hamiltonian 



p,q=i 



+ 3m(B pq ) i (flfl+fpfq) + 3m(A pq ) i ($f q +f p f}) . 

Note d\z{A) = SKe(A)*, Jm(A) = -3m{A)K SKe(B) = 
— 9 £ le(i?) t , and 3m(B) = — 3m(i?) t which is a consequence 
of A — A^ and B = —B f . We rewrite the Hamiltonian 
using Majorana operators such that 



XX ■ ■ II + YY • ■ II 



II • • XX + II • ■ YY 



-iH = — 



^2 -y\t{Ap P ) m 2 p-im 2 p + ^ V p< 



L P =i 



p,q=l;p>q 



where 

Vpq = - Ule{B pq ) (m 2 p-im 2 , - m 2q _im 2p ) 

- y\t{A P q) {m 2 p-im 2q + m 2q -im 2p ) 

- 3m{B P q) (■m 2p -im 2q _ 1 - m 2 pm 2q ) 

- 3m(A pq ) (m 2 p_ 1 m 2 q-i + m 2p m 2 q) . 

By applying the Jordan- Wigner transformation we ob- 
tain the Hamiltonian for the corresponding spin system 
(for better readability, the tensor-product symbol is omit- 
ted, e.g., IXY :=I®X® Y): 



a a 



j—l d—p 



• ,q=l;p>q 



on n spin-^ qubits with n > 2. Assume that one end qubit 
is individually locally controllable. The system algebra is 
given as the subalgebra so(2n + 1) which is irreducibly 
embedded in su(2 n ). 

Proof. We use the fermionic picture where the number of 
modes d equals the number of spins n. The generators are 
given by w\ = L(v\) with v\ = J2'p=i - m 2p-im 2p+2 + 
m 2p m 2p +i, L(mi), and L{m 2 ). Obviously, the element 
w 2 = L(v 2 ) with v 2 = m,\m 2 can also be generated. 
One can verify that exactly all Majorana operators of 
degree one or two can be obtained: One line of reason- 
ing uses Lemma [T4l together with the commutation rela- 
tions [L(m 2p -i), (bp)} = L(m 2p+2 ) and [L{m 2p ) , L(b p )] = 
—L(m 2p+ i) to show that all degree-one operators can be 
generated. This immediately gives all quadratic oper- 
ators as well, while operators of higher degree are not 
attainable. Therefore, the dimension of the system al- 
gebra is 2c? 2 + d. Note that the operators L(m 2p -im,2p) 
form a maximal abelian subalgebra which proves that 
the system algebra has rank d. One can show by explicit 
computations on the generators of the spin system that 
the system algebra is embedded irreducibly in su(2 d ). As 
the coupling graph of the spi n sy stem is connected, we 
conclude with Theorem 6 of [18| that the system alge- 
bra is simple. Listing all simple (and compact) Lie alge- 
bras with the correct dimension and rank, we obtain (a) 
So(2d+l) for d > 1, (b) sp(d) for d > 1, (c) su(2) = so(3) 
for d = 1, and (d) z§ for d — 6. As the system algebra 
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contains also all quadratic operators, it has a subalge- 
bra so (2d) which is of maximal rank. This rules out the 
cases (b) and (d) (see, e.g., p. 219 of [6(J)- And for d = 1 
the cases (a) and (c) are isomorphic. This completes the 
proof. □ 

Note that with our fermionic approach one can readily 
determine the dimension and rank of the system alge- 
bra. Likewise, we establish that all fermionic operators 
act irreducibly from which we can infer that the system 
algebra is simple. The rest of the proof follows by an 
exhaustive enumeration. — In more general terms, as in 
Theorem 34 and Corollary 35 of \Xf$, we obtain a con- 
nection between a spin system and a fictitious fermionic 
system: 

Corollary 45. Consider a fictitious fermionic system 
with d modes which consists of all linear and quadratic 
operators and whose generators can, e.g., be chosen as all 
Majorana operators of type L(m 2 p—i) combined with the 
Hamiltonian from Eq. (|21[) where the control functions 
.A.pq Q/nd Jl3'pq can be assumed to be real. This fictitious 
fermionic system and the spin system of Proposition \44\ 
with n = d spins can simulate each other. In particular, 
both can simulate a general quasifree fermionic system 
with d modes and system algebra so (2d) as presented in 
Proposition \13\ and Theorem \15\ 

2. A Spin System with System Algebra so(2?i + 2) 



Principle Remark: Now we have established a setting 
that allows for exploiting the powerful general results of 
[73| on the structure of orthogonal groups that provide 
a second avenue to Proposition 03] assuming we have al- 
ready established Proposition 0HJ Lemmata 3 and 4 of 
(73l | show that for k > 3 any subalgebra of so(k) with 
dimension (k — l)(k — 2)/2 is isomorphic to so(fc — 1); 
moreover so(fc — 1) is a maximal subalgebra of so(fc). 
Thus, by proving that the system algebra has dimension 
2d 2 + d with d > 1, it can be identified as the subal- 
gebra so(2d + 1) of so(2d + 2). We emphasize that this 
particular proof technique should be widely applicable in 
quantum systems theory. 

Relying on the proof of Proposition |46| and building 
on Theorem 32 as well as Corollary 34 of [18|], we obtain 
connections between a spin system, a quasifree fermionic 
system, and a fictitious fermionic system: 

Corollary 47. The following control systems all have the 
system algebra so(2fc + 2) and can simulate each other: 

(a) the spin system of Proposition with k spins, 

(b) the quasifree fermionic system with k + 1 modes as 
presented in Proposition \13\ and Theorem \15[ and 

(c) a fictitious fermionic system with k modes which con- 
tains all Majorana operators of degree one, two, 2k— 1, 
and 2k, and whose generating Hamiltonian can be chosen 
from Eq. (|2ip where the control functions A pq and B pq 
can be assumed to be real. 



Proposition 46 (see Proposition 28 in [18j). Consider 
a Heisenberg-XX. chain with the drift Hamiltonian 

H d = XX ■ • II + YY ■ ■ II + • • • + II ■ • XX + II ■ • YY 

on n spin-^ qubits with n > 2. Assume that each of the 
two end qubits is individually locally controllable. The 
system algebra is given as the subalgebra SO (2n + 2) which 
is irreducibly embedded in su(2" ). 

Proof. We switch to a fermionic picture where the num- 
ber of modes d equals the number of spins n. The gen- 
erators are wi = L(v\) with v% = J2pZi — m 2p~im 2 p+ 2 + 
m 2p m 2 p +1 , L(mi), L(m 2 ), L(m 2 d-i m 2p -im 2 p), 

and L(m 2 dY[p = i m 2p-i'm 2 p). One can verify by explicit 
computations that exactly all Majorana operators of de- 
gree one, two, 2d— 1, and 2d can be generated. Therefore, 
the dimension of the system algebra is 2d 2 + 3d+ 1. Note 
that the operators L(m 2p -im 2p ) together with the oper- 
ator i(Op=i m 2p-inrip) form a maximal abelian subalge- 
bra which proves that the system algebra has rank d + 1. 
One can show by explicit computations on the generators 
of the spin system that the system algebra is embedded 
irreducibly in su(2 d ). As the coupling graph of the s pin 
system is connected, we conclude with Theorem 6 of |18j 
that the system algebra is simple. The proof is completed 
by listing all simple (and compact) Lie algebras with the 
correct dimension and rank. We obtain (a) so (2d + 2) for 
d > 1 and (b) su(4) = so(6) for d = 2. □ 



Appendix C: Proof of Theorem 1171 

The cases of d G {2,3,4} can be verified directly. We 
build on Lemma [14] and obtain a basis of t\ consist- 
ing of L{a p ) with 1 < p < d as well as L(bp^) with 
bp := -m 2p _im 2p+2l + m 2p m 2p+2i -i and L(cp ) with 

cy := m 2p ^im 2p+2i -i + m 2p m 2p+2l where p, i > 1 and 
p + i < d. One can systematically enlarge the index 
(i) starting from the elements L(& p ^) = L(b p ) G ti and 
L(cp^) = (c p ) £ 61 and generate all L{b p ^) and L(c p ^) by 
combining the commutator relations [L(c p ), £(6^^)] = 
-L(b p l+1) ) and [L(c p+l ),L(b p l) )} = L{b p l+1) ) with the 
commutator relations [L(a p ), L(b p ^)] = — L(c p ^) and 

[L(a p ),L(cp )] = L(b p ). It is straightforward to check 
that no further elements are generated by commutators 
starting from the elements L{a p ), L(b p ^), and L(c p ^). 
We obtain that dim(li) = d + (d - l)d = d 2 . Fur- 
thermore, the elements L(a p ) form a maximal abelian 
subalgebra of ii and the rank of ti is d. It follows 
that £1 is a subalgebra of maximal rank in so (2d). We 
now show that the center of fii is one-dimensional and 
is generated by L(c) with c := X)p=i a p- Combining 

the commutator relations [L(a p +i ),L{b { p l) )]=L{c p l) ) and 
[L{a p +i), L(c p ')] = — L(b p •*) with the ones for L(a p ) men- 
tioned above, we conclude that [L(a p +a p +i), L(b p ■*)] = 
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[L(a p +a p+ i) 7 L(cp )] = 0. In addition, we obtain 

[L( aj ),L(b^)} = [L(a 3 ),L(c p l) )] = if p jt j jt p + i. 

It follows that [L(c),L(b$)] = [L(c), L{c p l) )} = and 
that L(c) commutes with all elements of ti. We rule 
out the existence of further elements in the center by 
explicitly computing the semisimple part s := [ti,ti] of 

ti. By applying [L(bj^), L(c p ^)}/2 = L(a p+i - a p ) com- 
bined with previously mentioned commutator relations, 
we can fix a basis of s consisting of the elements L(b p ^), 
L(c p ^), and L(a p — a p+ \) where 1 < p < d — 1 and 
1 < i < d — p. We proceed to prove that s is actu- 
ally simple by showing that s is not abelian (which ob- 
viously holds) and that any non-zero ideal i of s is equal 

to s. Starting from (ad 2 (L(a g ))) L(b p ^) = —L(b p ^) and 
(ad 2 (L(a q ))) L(c p ^) = -L(cp ) ) for q = p or q = p + i, 
we deduce that ad 2 (L(a q )) + ad 4 (L(a q )) = 0. Like- 
wise, y p ) := [ad 2 (L(a p -a p+i )) + ad 4 (L(a p -a p+l ))]/12 
annihilates all basis elements of s except for L(b p ^) 
and L(cp ) which are left invariant. Using the defini- 
tion xf := [&& 2 {L{b p l) )) - ad 2 (L(4 l) ))]/4 and verify- 
ing Xp L(a p ) = Xp L(a p+ i) = 0, we can infer that 
x p L(a q —a q+ j) = holds for all valid q and j. Fur- 
thermore, we have Xp L(b q ^) = Xp L(c q ^) = for all 
valid q and j unless when both q = p and j = i hold. We 
obtain x p l) L{b p l) ) = L{b p l) ) and x p l) L{c p l) ) = -L{c p l) ) 
in this exceptional case. As s is semisimple, i cannot 
be abelian and has to contain an element which is sup- 
ported on some L(b p ) or L(c p ). Relying on the ideal 
property [s,i] C i and the operators x p and y P \ we 
conclude that L{b p ^) 6 i or L{c p ^) <G i. Obviously, the 



conditions L(b p ) € i, L{c p ) G i, and L(a p —a p+ i) £ i are 
equivalent. By applying previously mentioned commuta- 
tor relations, we can verify that L(b q ^) <G i holds for all 
q < p and q + j > p + i. In particular, L(6^ _1 ) £ i. Us- 
ing the commutator relations [L(c p ), L(b p ')] = L(b p 1 ^) 
and [L(c p +i_i), L(b p ')] = —L(b p 1 ') where i > 1, we 

can reach the conclusion that L(b q ^) S i for all valid q 
and j. Thus, we have shown that i = 5 and s has to 
be simple. We summarize that ii has dimension d 2 , has 
rank d, and is a subalgebra of maximal rank in so (2d). 
In addition, it is a direct sum of a simple subalgebra 
and a one-dimensional abelian subalgebra. All maxi- 
mal subalgebras of maximal rank in so(2d) are given 
by (a) so(2fc) + so(2d - 2k) with 2 < k < d - 2, (b) 
so(2d — 2) + u(l), and (c) u(d) = su(d) + u(l) (see, e.g., 
p. 219 of [6(3]). All other subalgebras of maximal rank 
in so (2d) are proper subalgebras of the parent algebras 
given in the cases (a)-(c); moreover they have to refine 
the direct sum decomposition of their parent algebra (see, 
e.g., p. 205 of [60]). The case (a) and its proper subal- 
gebras are ruled out by the fact that their semisimple 
part would not be simple. The proper subalgebras of 
(b) can be analyzed by iterating the possibilities (a)-(c) 
with respect to so (2d — 2), and one concludes that ei- 
ther their center would not be one-dimensional or their 
semisimple part would not be simple. All the proper 
subalgebras resulting from case (c) have a dimension less 
than d 2 . We can have a dimension of d 2 for the algebra 
so(2d— 2) + u(l) of case (b) only if d = 4 when (b) and (c) 
coincide as so (6) = su(4). Thus (c) is the only possibility 
left, which finally completes the proof. 
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